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Chapter 1 

Introduction 



One of the principal goals of physics is to understand the constituents of matter 
and their interactions with each other. Experimental observations provide the 
basis for theoretical models, and these models must in turn be tested by more 
experiments. Theoretical physicists invariably introduce certain parameters to 
simplify these models. For example, the concept of electric charge is introduced 
to find the relative electric force on the object. An object with twice the charge 
will experience twice the force. Charge is a fundamental quantity and can be 
used with complete ignorance of the inner structure of an object. 

We arc interested here in another fundamental, yet less-understood, quan- 
tity: polarizability. The origin of this concept is in the study of atomic spectra. 
It was observed that placing atoms in a static electric field lowers their energy; 
the polarizability was introduced to help quantify the amount of this energy 
shift. Polarizability is important because it is a basic property of a particle, 
along with charge and mass, and is involved in the description of all electro- 
magnetic interactions. 

Part of the interest in this subject stems from the fact that as an experimen- 
tal quantity, the polarizability is a single number which represents the influence 
of many effects which cannot be directly measured. However, these effects can 
be calculated in various theoretical models and compared with the experimen- 
tal measurements. In the quark model, for instance, a proton is composed of 
three smaller particles known as quarks. One of these is negatively charged and 
the other two are positive. An electric field would exert a force on the positive 
charges in the opposite direction as the negative charge, causing the proton 
to deform slightly. The amount of this deformation depends on, among other 
things, how strongly the quarks interact with each other. The polarizability, 
on the other hand, contains information on the deformation, but no knowledge 
of the quark interactions is needed to measure it. It can therefore be used to 
independently check the various models which have been proposed to describe 
nuclear matter. 

The focus of this paper is on the polarizability of the neutron. To date, 
accurate measurements of the neutron polarizability have not been realized. 
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This is mainly because neutrons have no charge, and so the interaction with 
electromagnetic fields is weak. The central question to be considered is whether 
the neutron polarizabililty can be determined in Compton scattering from a 
dcutcron, which is a particle containing one proton and one neutron. This allows 
the use of a stable target. Unfortunately, the computation is more difficult for 
the deuteron because all of the proton interactions also need to be taken into 
account. As discussed in the following section, this is not a new idea and there 
have been previous calculations. However, advances on both the theoretical and 
experimental fronts invite the possibility of not only a more accurate calculation, 
but a more accurate measurement. The intent of this work is provide a complete 
calculation which would be useful to an experimentalist interested in measuring 
the neutron polarizability. 

Chapter 2 contains an introduction to polarizability and examines the pre- 
vious theoretical and experimental work which has been done in this area. De- 
scriptions of the methods used in the calculation are discussed in Chapter 3. 
Checks which were performed to ensure the correctness of the work are also 
detailed here. The actual results of the numerical calculation are presented 
in Chapter 4, along with the discussion of the question of neutron polarizabil- 
ity measurements from deuteron Compton scattering. The conclusion can be 
found in Chapter 5. Details of calculations omitted in the main text, as well as 
formulas for all of the relevant scattering amplitudes, are in the Appendices. 



Chapter 2 

Background and Motivation 



2.1 Definitions of Polarizability 

Polarizability, like charge, is one of the fundamental properties of a particle. 
Roughly speaking, it is a measure of how easy it is to deform the object in the 
presence of an electric or magnetic field. The symbol a is traditionally used to 
represent the electric polarizability, and /? is used for the magnetic polarizability. 

A more precise definition can be given in terms of induced dipole moments. 
The polarizability is the constant of proportionality between an external field 
and the average dipole moment that it induces. We can write 

((?) = aE, (2.1) 
(/?) = (2.2) 

where (d) ((/i)) is the average induced electric (magnetic) dipole moment. This 
induced moment in turn changes the potential energy by an amount 

AE = -ia£;2 _ 1/3^2 (•2 3^ 

This energy shift is known as the quadratic Stark effect. 

We can also formulate a matrix definition of polarizability. The Hamiltonian 
for a dipole interacting with an external field is given by 

ifint = -d,E, (2.4) 

where we have chosen E to be in the 2;-direction for simplicity. The energy 

shift can be calculated using perturbation theory. The first-order term is zero, 
since the dipole operator changes parity. Restricting ourselves to nucleons, the 
second order term is 

Em — Emi 
N' " " 



8 



CHAPTER 2. BACKGROUND AND MOTIVATION 



9 




(a) 



(b) 



Figure 2.1: Feynman diagrams for Compton scattering, uncrossed (a) and 
crossed (b) terms. 



where iV is a ground-state nucleon and N' is some excited state. An upper- 
case N is used throughout the text to denote a generic nucleon, and n and p 
represent the proton and neutron specifically. Comparing this expression with 
equation (|2.3|), we see that the polarizability can be written as 



N' 



Similarly, 



N' 



\{N\d,\N')f 
En' — En 



\{N\fi.\N')\' 
En' — En 



(2.6) 



(2.7) 



Thus, polarizabilities can be calculated by taking matrix elements of dipole 
operators between ground and excited nucleon states. 

We can make a rough order-of-magnitude estimate of by letting N' be 
an A^* state. Then, 



aN 



2e2|( A I z I A* )|2 



20 X lO^'^fm-' 



TTlN' 



rriN 



(2.8) 



where the square of the matrix element was estimated to be one-third of the 
mean nucleon radius squared, or about ^fm^. This is on the same order as 
the currently accepted value of aN « 12 x 10^'^fm^. Hereafter the unit of 
polarizability will be assumed to be 10~*fm'^ . 

The above expressions define the polarizability in terms of either of an en- 
ergy shift or an induced dipole moment. A third definition, which turns out 
to be the most practical, involves the scattering amplitude for elastic photon 
scattering, commonly known as Compton scattering. The Feynman diagrams 
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for this process are given in Figure 2.1. The initial photon has four- momentum 



{uji,ki) and polarization vector cx^, where Xi — ±1, while {u}f,kf) and eaj- de- 
scribe the final photon. The angle between ki and kf is labeled 9, and the initial 
and final photon energies are related by the kinematic relation 

- l-^^(?-cos0) • (2.9) 

TUN ^ ' 

The scattering amplitude (to order w^) can be written as ||] 

e2 

fiuj^uji) = {el ■ £xj + it^i+ujf)g{e,mN,KN;kf,el h,ix,) + 

aujfUJi{e*x^. • caJ + (3ujfuj^{kf x e^^) • {h x eA J + 

ujfUih{e, TUN, kn; kf,i*x^ , fc^, ca^)- (2-10) 

The first term in the amplitude, known as the Thomson amplitude, is the dom- 
inant term at low energies. It depends on only two parameters: the charge and 
mass. Since neither of these quantities assumes any internal structure for the 
nucleon, this is consistent with the fact that in scattering at low energies the 
nucleon must behave like a spinless point particle. 

The next order term, represented above by g, is a function of not only the 
charge and mass, but also the anomalous magnetic moment k. Therefore, to 
this order, the nucleon still acts like a point particle, but one with spin. Not 
until the second order does any internal structure become important. This 
piece of the amplitude is divided into two parts. The terms with the angular 
dependence (e^^ • caJ and {kf x e^^) . {ki x ixj are written separately, and all 
other combinations of the four photon unit vectors are lumped together in the 
function h. We can define the coefficients of these two terms to be the generalized 
polarizabilities a and p. This is the most general definition of polarizability. By 
defining a and f3 in this manner, they can be determined experimentally by 
examining the angular dependence of the cross-sections at higher energies. 

The Compton scattering definition involves one additional complication: the 
generalized polarizabilities a, (3 are not the same as the static polarizabilities 
a,fi. They are related by 

a = Q;-hAa. (2.11) 

Since the polarizability is defined as a coefficient in the scattering amplitude, 
there is no way to experimentally distinguish between different terms which have 
the same energy and angular dependence. The main contribution to a comes 
from the effects of the induced dipole moments, and all other contributions are 
collectively called Aa. 

The Aa term includes corrections which account for the finite size of the 
nucleon. The first term in the explicit expression for Aa is ^ : 

Ao^jii + 5M, ,2,12) 

6ra Am"^ 
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where Ve stands for the charge radius. Aa accounts for about 40% of ctp, while 
it is only about a 5% effect in the neutron. This can be traced to the fact that 
the neutron has zero charge. 

Corrections also exist for the generalized magnetic polarizability. The first 
term in A/3 is M 

e p 

AP^-^. (2.13) 

Additional corrections involve the nucleon form factors i^i(4r7i^), F2(4m^) and 
their derivatives 



2.2 Model Calculations of Polarizability 

While we are mainly interested in determining polarizabilities through exper- 
iment, it is useful to compare these values with those obtained in theoretical 
calculations. Many attempts have been made to obtain values for a and /3, 
but all of the results are model-dependent. We will examine some of these re- 
sults, including those obtained from dispersion relations and chiral perturbation 
theory. 

One of the first to obtain an estimate for the proton polarizabilities was 
Baldin jij. His calculation of 

4.0 <ap< 15.0 (2.14) 
is still consistent with recent calculations. He considered only the first excited 



state {N' — N + tt) in the matrix definition of Op (equation 2.6), and estimated 
the matrix elements using the meson photoproduction data available at the 
time. 

Several useful dispersion relations can be derived for the generalized po- 
larizabilities. Gell-Mann and Goldberger Q wrote down the once-subtracted 
dispersion relation for the forward scattering amplitude /(w): 



X ./X lmf(Lu')duj' 



Looking back at equation ( ^.10 ), we see that the actual form of real part of 



the forward scattering amplitude is 



e2 



Re/(u;) = + {a + l3)Lo\ (2.16) 

since the function g is purely imaginary, and contributions from h cancel (this 
can be seen from the explicit expression for h'm. ||l|). After comparing the above 
equations and then using the optical theorem, 

Im/(L^) = ^atot(^), (2.17) 
47r 
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where crtot(w) is the total cross-section for photoabsorption, we can write the 
dispersion relation Q 

CTtot {i^')duj' 



This relationship, the Baldin sum rule, constrains the possible values of a 
and /3. It is model- independent, and experimental data alone can be used to 
evaluate the integral, using some reasonable method to continue to infinity. The 
generally accepted results from the sum rule are 0] 

ap + (3p = 14.3 ±0.5, (2.19) 
an+(3n = 15.8 ±0.5. (2.20) 

However, a recent experiment predicts the following: 

ap + Pp = 13.69 ±0.14, (2.21) 
a„+(3„ = 14.44 ± 0.69. (2.22) 

The values for the proton are consistent with the older ones, but the new neutron 
values are slightly lower. 

A dispersion relation for a — f3 can also be derived |9[ ^ , but it contains a 
contribution that depends on the amplitudes for 77 — s- tttt, which are not well 
known. Nevertheless, this difference has been estimated to be 

ap-/3p « 3.2, (2.23) 
a„-Pn - 3.9. (2.24) 

Comparison with the sum rules above yield smaller values for the individual 
polarizabilities than the « 11 and a„ ~ 12 measured experimentally. 

Various quark models have also been used to determine polarizabilities. A 
simple nonrelativistic model with a harmonic oscillator potential yields reason- 
able values for f3 but predicts dip > q;„, so more sophisticated models are needed. 
For example, Werner and Weise jl^ calculated a and /3 using a valence quark 
core surrounded by a pion cloud to be 

aN ~ 7-9, (2.25) 
Pn ~ 2. (2.26) 

Again, these values are lower than recent experimental measurements. 

Perhaps the most promising method is using chiral perturbation theory to 
calculate the polarizabilities. Calculations with only the one-loop contribution 



produce |13 



aN = 10P^ = lf^—^ «13.6. (2.27) 

Including the next order as well as the effects of the A resonance gives [|l^, |l^ 

Op = 10.5 ±2.0 a„ = 13.4 ±1.5 
/5p = 3.5±3.6 /3„ = 7.8±3.6 
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If we compare these two calculations with the Baldin sum rule, we find agree- 
ment with the predictions for the proton. However, the earlier calculation gives 
better agreement with the neutron than the higher-order one, indicating a need 
for more work on this front. 

Both of these methods make predictions about the sources of the various 
contributions to the polarizability. For example, it can be shown that the most 
important contributions to a come from polarizing the pion cloud - about 50- 
70% ||l^. Resonances such as the A and the pion polarizability account for 
the remainder. The small value of /3 w 2 can be attributed to cancellations 
between the large paramagnetic contribution from the A resonance and the 
diamagnetic contribution from the pion cloud (A/5) |p^ . 



2.3 Proton Experiments for Determining Polar- 
izabilities 



We now examine some experimental results. Several experiments have been 
performed over the past few years with the aim of measuring nucleon polariz- 
abilities. Most of these have specifically targeted the proton, but some attempts 
have also been made at measuring the neutron polarizabilities. Since there have 
been so many more measurements of ctp and Pp, we begin with the proton. 
The differential cross-section for proton Compton scattering, which is pro- 



portional to the square of the scattering amplitude of equation (2.10), can be 



written to order lu^ as 



da da 



1 



— = — uju;' \-{ap + Pp)il + cos9f + -{ap-Pp){l- cos ey 

dil ail Born nip \u} J [2 2 

(2.28) 

where a is the fine structure constant, and the Born term, first written down 
by Powell 0, is 



da 

rfil Born 



1 

2 \mp 

2LJiL0f 

•ip n 

.3 ^^^f 



2 / \ 2 



(1 — cos( 



— ^ - 

bJf LUi 

"T i^p 5~ 



4(1 - cos( 



1 



(1 — cos ( 



kI-^ [2(1 -cos^) +sin2 



^UJiUJf 



1 



(1-H-sin^ 



Note that the polarizability terms in equation ( 2.28|) arise from interference 
between the polarizability amplitude of equation ( ^.1C| ) and the Thomson am- 
plitude. 

This formula illustrates the sensitivity of the cross-section to a.p and (3p at 
different angles. The cross-section at forward angles is most sensitive to Up + ^p, 
while at backward angles only oip — (3p can be measured. At 90° the /3p terms 
drop out completely. 
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Another point to consider is that there is an optimal range of energies for 
which a measurement of the proton polarizabihty can be done. The expansion 
to order is only valid to about 100 MeV ||l^. On the other hand, since 
the polarizability terms are of this order, too low an energy means that the 
cross-section is too insensitive to cip and f3p. Figure 2.2 shows the dependence 
of the polarizability terms on energy for two different angles. Looking at this 
graph, it seems that the optimal energy range for an experiment would be 70-100 
MeV in order to balance sensitivity with expansion validity. The cross-section 
also depends more heavily on the polarizabilities at forward angles; however, 
these measurements alone would only test the Baldin sum rule and not yield 
individual values for a and /3. 

A summary of the data from a few recent proton Compton scattering ex- 
periments is given in Table below. The Moscow experiment ||l^ is included 
because it was believed to be the most accurate of earlier attempts to measure 
the polarizabilities, despite the low energy resolution of the photon detector. 
Their value for cip is consistent with theoretical predictions, but ap + j3p falls far 
below the accepted value from the sum rule, and Pp is negative. This formed 
part of the motivation for more recent experiments . The values listed in the 
table have actually been generated by a recent re-analysis of the original data 

[io|. 

The Illinois group performed one of these latest experiments. Measure- 
ments were taken at both a forward and a backward angle, and the energies used 
fall well within the range of validity of the low-energy expansion. However, this 
was done at the expense at sensitivity to the polarizabilities, which explains 
the high uncertainties. An improved 1995 experiment by the same group pO| 
used higher photon energies, but now perhaps too high to permit use of the 
low-energy expansion. 

A 1992 Mainz experiment |Q is also included in the table. Since data was 
only taken at the backward angle, only the difference OLp — (3p could be extracted. 
Individual values for oip and (3p were obtained from the sum rule. Also, the 
energy range used is possibly beyond the limit of validity of the low-energy 
expansion. 

More extensive reviews of these and other proton Compton scattering exper- 
iments can be found in the literature |jl|, |l^. However, the point to be made is 
that, despite some shortcomings, measurements for the polarizabilities in these 
experiments are in reasonable agreement with each other and with theoretical 
predictions. This is not the case for the neutron. 



2.4 Neutron Experiments for Determining Po- 
larizabilities 

The neutron polarizabilities have not yet been measured with reasonable accu- 
racy. This is mainly because a neutron Compton scattering experiment cannot 
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Figure 2.2: Energy dependence of a and /3 terms 
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Table 2.1: Experimental results for ap and Pp 



Wi(MeV) 6* (degrees) 



Moscow(1975) 
Illinois (1991) 
Mainz(1992) 



70-110 
32-72 
98, 132 



90,150 
60,135 
180 



Illinois/SAL(1995) 70-148 90,135 



result 

ctp + jSp — 5.8 ± 3.3 
ap- Pp = 17.8 ± 2.0 
ap = 10.9 ± 2.2 ± 1.3 

Pp = 3.3 T 2.2 =F 1-3 



R — 7 nq+2-49+2.14 
IJp — ' •^J'J_2 37_2.05 
— 1 n co-(-l-25+1.07 

o ro + 1.25+1.07 
"'■^°-1.19-1.03 

ap = 12.5 ±0.6 ±0.7 ±0.5 
(3p = 1.7 T 0.6 T 0.7 T 0.5 



JLp 

ap 
Pp 



be directly performed. Not only is a neutron target unavailable, but the cross- 
section for this process would be extremely small. Therefore, other methods of 
measuring a„ and /3„ have been developed. The two that have produced results 
thus far are Coulomb scattering and deuteron photodisintegration. 

Coulomb scattering of a low-energy neutron off a heavy nucleus is a reason- 
able means of measuring the electric polarizability because of the large dipole 
moment induced by the Coulomb field. The interaction energy is given by 

1 Z^e^ 

V = --a,,—^, (2.29) 

where Ze is the charge of the heavy nucleus. Note that only q;„ (rather than an) 
is involved here, and the Aa contribution, although small, must be accounted 
for. 

This energy shift turns out to have only about a 1% effect on the cross-section 
[ p3| . The low sensitivity to a„, combined with the uncertainty in calculating 
other small contributions to the cross-section, makes an accurate determination 
of the polarizability in this type of experiment very difficult. Two separate 1988 
experiments, both using ^ospj^ ^j^g ig^^get nucleus, quote large uncertainties 




an = 12.0 ±10.0 (Vienna/OakRidge), (2.30) 
an = 8.0 ±10.0 (Munich). (2.31) 

A more recent experiment by Schmiedmayer et al. | |2^ , also using a lead 
target, produced the following result: 

a„ = 12.0 ±1.5 ±2.0. (2.32) 

This is the value normally quoted for the neutron electric polarizability. How- 
ever, the accuracy of this result has been questioned |2^, |2|] . It has been claimed 
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that the uncertainties were underestimated, and that the best estimate for a„ 
is reaUy ~ 7 - 19 

The other method used to measure a„ is quasi-free Compton scattering by 
the neutron bound in the deuteron. It has been calculated that the cross-section 
for deuteron photodisintegration is sensitive to the neutron polarizability when 
most of the photon's momentum is transferred to the neutron, and the proton 
remains in the target pg[ |. Such an experiment was carried out by Rose et al. 
in 1990, and yielded the measurement 

an = n.7ttf,r (2-33) 

Again, the uncertainties are large. However, there have been recent arguments 
that this method should be revisited as it is capable of producing more accurate 
resuhs |l|. 

It is clear that better experimental measurements of the neutron polarizabil- 
ity are needed in order to keep up with recent theoretical developments. This is 
why we have chosen to investigate deuteron Compton scattering as a means of 
obtaining additional information about q:„. This idea seems to have been orig- 
inally suggested by Baldin who calculated the cross-section in the impulse 
approximation. A more extensive calculation was undertaken by Weyrauch 
[HI but certain deficiencies in his method to be discussed later, as well as 
a lack of emphasis on the neutron polarizability, suggest that this question can 
be reexamined. 

A deuteron Compton scattering experiment was performed in 1994 p^ , 
using 49 and 69 MeV photons, and experiments are in progress at Saskatoon 
[p5| and Lund [p6[ . This has spurred some recent activity in this area. Articles 



by Arenhovel ||3j and L'vov and Levchuk 39 quote new theoretical results 
for Compton scattering. A calculation using effective field theories has also been 
published Their results will be compared with the present calculation in 
Chapter 5. 



Chapter 3 



Deuteron Compton 
Scatttering Calculations 



3.1 Feynman Diagrams 

The Feynman diagrams for several of the processes contributing to deuteron 



Compton scattering are shown in Figure 3.1. The ovals at either end represent 
the deuteron wavefunctions, while the solid lines are the individual nucleons 
with which the photons (wavy lines) interact. Not all possible combinations 
are shown; interactions can occur on either nucleon with either an incoming or 
outgoing photon. 



Figure 3.1(a) is known as the seagull diagram. This is the interaction that 
arises from the term in the Hamiltonian that is proportional to A^. It has 
no overall energy dependence, and is an important term at all energies. All 
other one-body interactions, which are at least of order w^, are represented by 



Figure 3.1(b). These include the terms which depend on the polarizabilities a 



and /3, as well as some relativistic corrections. 

Dispersive diagrams without meson exchange are depicted in Figure 3.1(c)- 
(d). At lowest order, each vertex can be either an electric or magnetic dipole 
(El or Ml) interaction, which originate in the p ■ A and jl ■ B terms in the 
Hamiltonian, respectively. These two diagrams account for the bulk of the work 
in the deuteron Compton scattering calculation. 

Three of the meson-exchange effects which are included in this calculation 



are shown in Figure 3.1(e)-(g). They are either second- or third-order terms, 



and are rather small at the energies considered here. The term of Figure 3.1(f) 
is called the "vertex correction" . 
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Figure 3.1: Feynman diagrams contributing to deuteron Compton scattering. 
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3.2 Overview of Calculation 

We will begin with an overview of the methods used in the calculation. The dif- 
ferential cross-section will be calculated non-relativistically using Fermi's golden 
rule, which can be written schematically as 

27r 1 2 

dafi = —6{AE) X . . — x (number of final states) x |7/j| . (3.1) 
ii imtiial nux 

The subscript fi indicates a dependence on the final and initial energies and 
polarizations. We use the convention where the initial photon flux is normalized 
to be 1 per arbitrary volume V. A detailed derivation of the phase space factors 
for this calculation can be found in Appendix D. We have set c = 1. 
The transition matrix Tfi is defined by 

Tfi = {df,Pf,^f\H"-'\di,Pi,7i)+ 
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E 

C,Pc 



df,Pf,jf I ff-t \C,Pc){C,Pc\ ff'"* I d„P„7, 



£;rf, + i^V2md + huj, - Ec 



P2 /2mrf 



IS 



(3.2) 



The notation is meant to represent all quantum numbers needed to de- 
scribe the initial (final) deuteron state, except for the center-of-mass momentum 
P, which is written separately so that it can be integrated out and absorbed 
into a momentum-conserving delta function. Similarly, C is the intermediate np 
(or possibly deuteron) state, and denotes all sums and integrals which 

are needed to describe this state. It is also possible for the intermediate state 
to contain photons and/or mesons, and these can be inserted into the above 
expression as needed. 

In order to calculate the transition matrix, we need to specify a Hamil- 
tonian. The basic non-rclativistic Hamiltonian, including polarizability and 
meson-exchange terms, is 



E 



^A\x,)-^A{x,)-pr 




e(l + Kj) 
2m -i 



a, ■ (V, X A{x,)) - (3.3) 



d^xe^A'^ix) [^+(f)(/)_(f) + (/)_(f)0+(x)] 



where 



^V. ) = ^ E 

fc,A=±l 

,(t) 



27rn, 



ik-Xj 



(3.4) 



The operators aj^ ^ destroy (create) a photon with momentum k and polarization 
A. The explicit time dependence e*"* has been suppressed since it only gives rise 



to the energy-conserving delta function in Fermi's golden rule (equation 3.1). 
Similarly, 



(j)±{Xj) 




(3.5) 



Here, a± destroys a ■k±. Finally, the tildes in the Hamiltonian indicate vectors 
in isospin space. 

The first term, as discussed before, is known as the Thomson term, and is 
the interaction represented by the seagull diagram of Figure 3.1(a). The next 
two terms are the electric and magnetic dipole terms. These two expressions 
will not be actually substituted into the transition matrix; instead, the more 
general observed vector current formalism, to be described later, will be used to 
describe the two-body terms. However, an explicit decomposition of the vector 
current expressions yields these dipole terms. 

The polarizability terms are next, and then the last three terms are the 
meson-exchange effects needed to calculate the diagrams in Figures 3.1(e),(f), 
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and (g), respectively. Relativistic corrections to this Hamiltonian will be dis- 



cussed in section 3.5 



3.3 Seagull Term 

The simplest diagram to evaluate is the seagull diagram, so it will be the first 
one we will consider. The basic procedure is to calculate the transition matrix 
using the interaction Hamiltonian 

i?^^ - ^A^ix-,), (3.6) 
Inserting this into the transition matrix gives 

^Z^'' = 9;;rv^^^ '^f^^f I • el e-^^r^^ \ d.,P, ). (3.7) 

^ ZlTlp V ^UJiUJf ' 

The next step is to eliminate the individual nucleon variables in favor of the 
deuteron variables r and R, where 

T = Xp Xyi^ (3.8) 

R = \{xp + x^). (3.9) 

We can now insert a complete set of and R' states into the matrix element. 
Defining q = kf — ki, the transition matrix becomes 

= /rfVd^i?'-:^^(^A, -A.) X 

^ J rUpV ^JUJ^ ^ 

( d/, F/ I e-t9-"*'e-*«-«~' I ?,R! )( ?^ R! \ d,, P, ). (3.10) 

The integral over R' can be performed by assuming that the center-of-mass 
wavefunctions for the deuteron are plane waves. This yields a momentum- 
conserving delta function. We remove this delta function for convenience by 
defining the scattering amplitude Mf^ by 



. 1 2^ ^W + ^,-P.-t.) 

■> V y^UJiUjJ V ■' 



Thus, 



= E / (^A, • eAj( I'SfJfMfVf I e-i'^-^ I r )( f\lSa^M,n ). 

(3.12) 
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Figure 3.2: Deuteron radial wavefunctions uo{r) and U2{r) 



Here we have explicitly stated what is meant by the label "di' , as the radial 
and angular parts have been separated. Isospin is not included as a quantum 
number, but is regulated by the generalized Pauli exclusion principle: 



(-1) 



L+S+T 



= -1. 



(3.13) 



Also, the orbital quantum number / is not written with an i or / label because 
the S and D states of the deuteron are summed at this stage. The spin and 
total angular momenta can only take on the values S = J = 1 for the deuteron. 

We now turn to a discussion of the deuteron wavefunction. The form that 
will be used here is 



OllAf ' 



U2{r) 



21 IM' 



(3.14) 



Instead of labeling the radial wavefunctions by the conventional u(r) and w(r), 
we use the notation ua{r) and U2{r) to create an explicit index which can be 
summed over. We use the wavefunctions published by Machleidt El[ (the "B" 



wavefunctions, derived from the Bonn potential). These are shown in Figure 3.2 
The angular wavefunctions are of the form ( f | LSJM ); an algebraic repre- 
sentation will not be needed. 



The easiest way to evaluate equation (3.12) is to expand the exponential into 
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partial waves using 



■ik-r 



J2 E ^^^''YLM{f)YEMik)j.{kr), 



(3.15) 



L=0 M=-L 



where jhikr) is the spherical Bessel function of order L. Inserting this along 
with the wavefunctions into equation ( |3.12 ) gives 



-in 



ex 



)E 

ll'LM 



u,{r)uAr){l'nMf \Ylm \ lllAU 



(3.16) 



This radial integral can easily be evaluated numerically; however, the matrix 
element must be simplified further. We go through this for the seagull term; 
derivations for other terms are in the Appendices. For the seagull term, the 
operator is a spherical harmonic, but other terms in the transition matrix will 
involve more complicated combinations of spin and angular momentum opera- 
tors. 

The first step is to form tensor products (each with a definite total angular 
momentum, which may be summed over) out of the individual operators. In 
this case, the operator has total angular momentum L. We can then use the 
Wigner-Eckart theorem (equation |K.2| ) to create reduced matrix elements which 
are independent of the azimuthal quantum number M. The operators in the 
reduced matrix elements can be uncoupled and evaluated using standard rela- 
tions (equations K.3-K.4), which can be found in a reference such as Rotenberg 
[ p2| . Further details of this for the seagull term can be found in Appendix A. 
The net result is: 



I'llMf I Ylm I lllM, ) = (-l)i-*^^+''3 



(2Z' + 1)(2L+1)(2/ + 1) 



47r 



i L r 





1 

-Mf 



L 1 

M M, 



1 I' 
I 1 



(3.17) 



Definitions of the 3j and 6j symbols can also be found in Rotenberg. One 

/ l„ 

important property to note here is that for the 3j symbol ( ^ ^ ^ 

I2, and Z3 satisfy the triangle inequality, and mi + m2 + m3 = 0. This 
means that, for the seagull term, M = Mf — Mi. Therefore, the sum over M in 



equation ( 3.16 ) can be removed and equation ( 3.17 ) can be inserted to produce 
the final result: 



m„ 



-(el,-eAj(-l)^-''^ X 



(3.18) 



1 

-Mf Mf 

-Mf Mf - M, M, ) ^2tAf/ -M. (9) (4°' + ^1" - ^) 
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where the notation 



has been used. 



Il^ / druAr)jU^)u,{r) (3.19) 



The largest term in equation ( 3.18 ) is the one for which I ^ V — L = 
0. This is because the amphtude of the radial wavefunction UQ{r) is several 
times larger than the amplitude for U2{r), as can be seen in Figure 3.2. Using 

I Ir n jir I = (—1)1+*^/ /V3, we can write the dominant part of the 

seagull term as 

=2 



Al^^ (dominant) = — (et • e^, )C<5A/,,Af. • (3.20) 
nip * ' 



This is the Thomson amplitude for proton Compton scattering reduced by the 
factor /go. 

The terms that depend on the nuclcon polarizabilities are very similar to the 
seagull term. Their transition matrices are calculated in detail in Appendices 
B and C. 

To find the differential cross-section, we insert the scattering amplitude into 
equation (D.12): 

^^ = (^\^ El I 12 



,o/ = - — l-^/.r, (3-21) 

ail J \ bJi J md 

where Ef = uji -\- md ojf, and is the deuteron mass. The spin-averaged 
cross-section is computed by summing over the final and averaging over the 
initial polarizations. 



3.4 Dispersive Terms 

Most of the remaining terms to be calculated are known as the dispersive terms, 
which are second-order in the interaction Hamiltonian. These terms, like the 
seagull term, are of order a, the fine structure constant. This is the small 
parameter in which we perform our perturbative expansion. 
The transition matrix has the form 

^ _ ^ \ {df,Pf,if\ g'"* \C,Pc){C,Pc\ g'"* I rf., P., 7. ) , 

^\ huJi + Ed,-Ec-P^/2md + ie + ^''■^^> 

C.Pc ^ 

{df,Pf,-ff I I C,Pc,j,,if ){C,Pcnulf I g'"* I d,,P,,-f,) \ 
-hujf + Ed,~Ec-Pl;/2md + ie J' 



Instead of inserting the Hamiltonian directly into equation ( 3.22| ), we make 



the substitution iJ'"' = — / J(^) • A{£^)d'^(^. To avoid confusion, ^ is used as a 
dummy variable, while x is a nucleon variable and r is a deuteron variable. 
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In addition, A is expanded into multipoles using 

L 



^Ae*«^=E E ^MA(0,-^,-^)^M^g^x (3.23) 



h^^^ + J-i^OyLMii) - iC^6.K)i"LA/(0 - ALyLM(|)j.(^0 

d'^\{0, — — (^) is the Wigner-d function (we use the convention of j43|) , 
which is the overlap of the state | LX ), rotated by the Euler angles (0, 
with the state \ LM ). This expression is derived in Appendix E. 

To simplify the notation, we will define a functions <i>i and <&/ by 



' eA.e*-« = V$,(0 + ---, (3.24) 




V$;(0 + ---. (3.25) 



The gradient terms above correspond to the first term in the brackets of equa- 
tion ( 3.23|) . In the low-energy limit, ^ ■ e\., so we expect e<I> to be 



responsible for the electric dipole interaction. 

Since A contains both a $ and a non-$ term, each i/'"* is the sum of 
two terms; in addition, there are four different possibilities for each term of 



equation (3.22). The largest term, which includes an El interaction at both Nj 
vertices, is the one for which both matrix elements contain It is calculated 
in detail in Appendix F, and certain aspects of it will also be discussed in this 
section. All terms in which exactly one of the matrix elements contain $ are 
calculated and tabulated in Appendix G. These include all magnetic and spin- 
induced interactions occurring at exactly one vertex. Most of the terms that do 
not contain $ are very small; only the results for the largest of these are listed 
at the end of Appendix G. 

As stated above, the largest second-order terms are the ones in which — J J(^)- 



V<I'(^)(i^^ is substituted for all four i?™*'s in equation ( 3.22 ). We want to avoid 
writing down an explicit expression for J, since this can only be done approxi- 
mately. Instead, we integrate by parts and use current conservation 

V-J(0 - -J [^''^(^^)] ' (3.26) 

PiO = Y.^M-^,)- (3-27) 

j 

The integral over ^ can then be performed inside the matrix element to give 

H'"^' = -J J{i) ■ V$(Od^f 

= -i[H,{e/h)^Xp)l (3.28) 
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where H is now the full Hamiltonian. If we examine only the piece of this 
commutator containing the proton kinetic energy, we see that 



2m,; 



(3.29) 



which follows from the canonical commutator [p, F{x)] — jVF{x) and the def- 
inition of $. This term is identical to one of the terms of the Hamiltonian in 



equation (3.3). Therefore, the $ commutators contain the operators responsible 
for electric dipole (El) transitions. 

We switch to the center-of-mass variables p= {pp — p„)/2 and P = pp + pn 
when we evaluate the commutator with the full Hamiltonian. Defining the 

2 

"internal" Hamiltonian H^^p = — h V, we get for the $i term 



if'"* = P-ex-e 



.iki-r/2^iki-B. _ 



^[i^"P,(e/n)$,(f/2)] 



(3.30) 



The term with the P operator will be called the "recoil correction"; its 
matrix elements are calculated in Appendix H. It is a small effect at the energies 
of interest (exactly how small will be discussed in section 3.3). The effects of 
all other recoil operators which occur in higher-order terms will be assumed to 
be negligible. 

The substitution Xp = r/2 + R must generate an exponential of the form 
gtP R^ which will combine will other similar exponentials to create the momentum- 
conserving delta function. This is the delta function which is removed in the 
matrix element M.fi that was defined in equation (3.11). The net result of ig- 



noring the recoil operators, therefore, is to replace Xp by r/2 and switch from 



The interesting physics is contained in the second term of equation ( 3.30| ). 
Plugging this into equation (3.22), we get 



M 



ft 




df I [H^P,^f]\C){C\[H-P,^,]\. 
huj., + Ed, -Ec- PSf^ma + ie 

df I [g"p,<l.] 

huJi + Ecu - Ec- Pll^mi ^ie \' 



(3.31) 



The dimensionless operators l>i = (e/?i)<i>i(r/2) and $/ = {e/K)'^ f{r/2) have 
been introduced, and the switch from 7^^ from M.fi has already been made. 
We have also labeled this term with the superscript "a" . 

Since | dij ) and | C ) are eigenstates of H^p , the commutators can be 
expanded and some cancellations with the denominator can be made. This 



generates four terms (after adding the crossed term of equation 3.22 



M}1 



(3.32) 
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where 



2md 



TlLOf + 



TlLOf + 



flLOj 



-| 2 

E 



df\<^f\C){C\<^.\d. 



hujj 



2mj 



Ed. -Ec+ie 



(3.33) 



2md 



2mj 



Edi - Ec+ie 



2md 



df I l>/4, I 



(3.35) 
(3.36) 



The first two terms correspond to the Feynman diagram Figure |3.l| (b) and 
(c), and includes the case where both vertices are El interactions. The third 
term is a small correction, calculated in Appendix F. The final term is respon- 
sible for compliance to the low-energy theorems which result from demanding 
gauge invariance. This is discussed at length in Section 3.6. 

The calculation of these terms, while more complex than the one-body terms, 
is straightforward except for the additional complications of the intermediate 
states and energy denominators. We use Green's functions to handle this prob- 



lem. Defining Eq = huji 
(neglecting overall constants) 



Ed 



the M'fl term can be rewritten as 



Eo 



■ le 



C){C\^,\d, 



(3.37) 



since we know that Ec is an eigenvalue of the -ff^c operator. The dependence 
of the Hamiltonian on the orbital angular momentum is now explicitly shown 
in order to remove the temptation to collapse the complete set of C states. 
However, if we split these states into angular and radial parts, the radial part can 
be removed, leaving a sum over the angular quantum numbers (LcScJcMc), 
which we denote collectively by C. Therefore, we split all parts of the transition 
matrix into radial and angular parts. The initial and final deuteron states can be 
separated according to the form of the wavefunction, equation (3.14). Looking at 
the algebraic representation of the <i> operator in the appendix, equations (F.2) 
and ( |F.3| ), we see that it also contains distinct angular and radial pieces, which 
we call O and J, respectively (suppressing the sums and related indices here 
for simplicity). Lastly, we insert two complete set of radial states to be able to 
write a meaningful expression: 

■^/' = EE E (^'ll^'^/|0/|(?)((7|0. MllAf,)x 

^oo pOO 

rr'ui,{r')Jj{r')G^^{r,r'-Eo).h{r)ui{r)drdr'. (3.38) 



CHAPTER 3. DEUTERON COMPTON SCATTTERING CALCULATIONS 29 



The Green's function is defined by 

G,^{r,r';Eo)^{r' \- | r ). (3.39) 

By applying the operator Eq — to both sides of this equation, the Green's 
function can also be seen to satisfy the following differential equation: 



h'ld' h'LcjLc + 1) ^. , , 



G,^{r,r';Eo) = 



S{r' - r) 



(3.40) 

For Vp(r), we use the Reid93 np potential as pubhshed in Q. 

All of the second-order terms can be reduced to the form of equation ( |3.38| ), 
so we need to be able to evaluate this double integral. Computing this integral 
directly would be a formidable task. Therefore, we evaluate it in two steps. 
First, we compute the function Xy-*^(?') defined by 

xf{r)^ Jr'dr'ui,ir')Jfir')G,^ir,r';Eo), (3.41) 

and then evaluate the single integral /"j*^: 



= / rdrui{r).Mr)x'f^{r). (3.42) 



The more difhcult step, of course, is the determination of (r). If we apply 
the operator Eq — to both sides of equation ( [3.41 ), we get 



c 



= ^iH^lEM, (3.43) 



This is an ordinary, inhomogeneous, second-order differential equation. We 
use the Numerov method |^ to obtain solutions for both the homogeneous 
and inhomogeneous solutions. However, the correct linear combination of these 
solutions still needs to be determined. The boundary condition is that x(^) 
must be an outgoing spherical wave at large distances. Assuming that J(r) — > 
as r goes to infinity, a condition always satisfied here, the boundary condition 
can be written as 

lim x(r) =4'^ (Qr), (3.44) 
where h'^l^ (Qf) is the spherical Hankel function defined by 



"Lc 

and where 



h^LliQr)^JLciQr)+^nLciQr), (3.45) 



= (3.46) 
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Q can be either real or imaginary. Therefore, to solve our problem, we must 
find a constant A such that 

lim {x'"'^°-(r) + Ax^^^'^lr)} = h^j^liQr). (3.47) 



Since the asymptotic solution of the differential equation ( |3.43 ) must be a 
linear combination of h'^i^^{Qr) and jLciQ''^)^ we can write the solutions for all 
r as 

Xir) = CLair)[jLciQr)+tLcir)h^LliQr)l (3-48) 

where the functions CLci'"') ^^^id tL^{r) approach the constants Clc t^^ as 
r approaches infinity. The choice of A which satisfies the boundary condition is 
then 

^inhom 

A --7^- (3-49) 

Our task is thus reduced to determining these coefficients. This must be 
done in the region where Cip(r) and t^ci''') are both constant, so that their 
derivatives vanish. We numerically evaluate x ^'^d its derivative, and then use 
the logarithmic derivative D = Inx(r) to solve for i^^, since 

, _ DjLc{Qr)-£jLc{Qr) 



After calculating t™^°™ and t^™, equation (3.48) can be used to determine 
^mhom g^j-^j Cjj°™. The value of A which satisfies the boundary condition can 
then be found. 



We have sketched the method used to calculate the integral in equation (3.41). 
In terms of numerical computation, this integral is the most demanding part 
of the calculation, so we would like to make a quick check on this algorithm. 
Since the Green's function is symmetric in r and r', the entire integral must 
be invariant under this interchange. Therefore, if / ~ V , making the switch 
Ji <-> Jf cannot change the integral, even if Ji ^ Jf. However, if it is true 
that Ji ^ Jf, the intermediate functions x('') created by the differential equa- 
tion routine will not be the same. There are several pairs of interference terms 
which have Ji ^ Jf and Ji Jf, so that comparing the values of the integrals 
for these terms would be a good check of the algorithm described above. This 
check has been performed on several pairs of terms, and the integrals are found 
to agree to 6 significant figures. 



3.5 Relativistic Corrections 

We now discuss relativistic corrections to our non-relativistic formulation. The 
major correction here is the spin-orbit effect. Corrections arising from boosting 
the final deuteron wavefunction into a moving frame are assumed to be small 
and therefore are neglected. 
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The origin of the relativistic spin-orbit effect can be understood by looking 
at the magnetic term in the Hamiltonian of equation (|3.3|): 



E 



2m i 



a, ■ B{xj). 



(3.51) 



The relativistic correction to the magnetic field for an object moving with ve- 
locity V is 

B^B-vxE, (3.52) 
This produces a relativistic term in the Hamiltonian: 



H 



RC 



E 

j=n,p 



e(l + 2Kj 



Ami 



■^3 ■ («j X 4) 



(3.53) 



The effects of Thomas precession have also been included. Since v ~ p/m and 
E = W{r) ~ r for a central potential, the dot product above goes as (? • L, 
where L is the orbital angular momentum. 
If we make the gauge invariant substitution 



and use 



H^'~^ becomes 



p — > p— eA, 



H 



RC 



E 



e(l + 2Kj) 



im'j 



(3.54) 
(3.55) 

(3.56) 



The leading relativistic correction to the scattering amplitude can easily be 
calculated from this Hamiltonian since it is a contact term, like the seagull term. 
Only the proton gives a correction because of the factor Cj . The result is 

- 2 

Mff - ( d/ I ■ {el^. X ex,){hLUf + hLu,){2Kp + l)e-'«-^"/2 \d, ). (3.57) 



Our calculation, however, does not use the Hamiltonian of equation (3.3) ; 
instead, the substitution H — — J J(^) • A(^)d'^^ was made and charge density 
and current operators were defined. To be consistent with this, we make a 
correction to the charge density: 



(3.58) 



where 

= eh 



^^^V^<5(e - r/2) ■ {ap x p) - ^V^Si^+ f/2) ■ (a„ x p) 
Ami mi 



(3.59) 
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This is derived from the non-relativistic reduction of the Dirac operators in 
the relativistic formulation. The corrections to the scattering ampHtude 
from this new operator are calculated in Appendix I. The leading-order term is 
identical to equation (3.57); higher-order terms, including one proportional to 
Kn, are neglected. 



3.6 Checking the Calculation 

There are checks which can be performed to help ensure that certain parts 
of the calculation are correct. One of these involves the use of the optical 
theorem to calculate the total cross-section for deuteron photodisintegration, the 
process in which a photon and a deuteron interact to form an unbound np state. 
These cross-sections have been more extensively studied, both theoretically and 
experimentally, so that the comparison of the value of the photodisintegration 
cross-section extracted from this larger calculation with the known values would 
be a good check. 

The differential cross-section for deuteron Compton scattering is given by 
(see equation D.12): 



(3.60) 



Combining the definitions of the transition matrix 7^^ (equation 3.1), the rela- 
tionship between 7^^ and -Mf^ (equation 3.11), and the mathematical identity 



1 



- le 



— h in6(x), 

X 



(3.61) 



we get a relation for the imaginary part of the forward scattering amplitude: 



lmMf,{e = 0) 



C I H" 



(3.62) 



where w is now the initial or final photon energy. To find the total photodisin- 
tegration cross-section dtot, we can use Fermi's golden rule. The result is 



ftot 



27r 1 ^ 

TiTf^ 



5{^E) I ( C I i?'"* I d ' 



(3.63) 



where C (as above) is an unbound np state. This can be combined with equa- 
tion ( 3.62 ) to yield 

47r 

atot = — ImX.,(e-0). (3.64) 
uo ■' 

This is the optical theorem. Since all of the jjS have already been calculated, 
determining the photodisintegration cross-sections is basically a matter of sum- 
ming the imaginary parts of the forward scattering amplitudes of all relevant 
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terms, and over all possible polarizations (of which there are only six, since the 
initial and final states are identical here). 

The only diagram that needs to be considered here is the uncrossed second- 



order term of Figure 3.1. However, all combinations of interactions at the two 
vertices must be included. In a 1964 paper, Partovi Q calculated the contri- 
butions from each of these interactions at 3 different energies. The dominant 
term is the one in which there is an El interaction at each vertex. He called 
this cross-section "approximation A" , and added the other terms one by one 
in successive approximations, until every term was included ("approximation 
I"). These results are reproduced in Table 3T, along with the corresponding 
cross-sections extracted from this calculation. There is reasonable agreement at 
all energies; any differences can be attributed to improvements in the wavefunc- 
tions and potentials (the short-range parts, in particular) since the time that 
Partovi's work was published. 

The meanings of the approximations in Table |3.l| : 

Approximation A: El only. 

Approximation B: A -|- singlet Ml. 

- E2. 

- triplet Ml. 

- spin induced triplet Ml. 
spin induced triplet M2. 

- spin induced triplet El. 
h retardation corrections to El. 

all other terms. 

We also compare these photodisintegration cross-sections with experimental 
data up to 100 MeV. This is shown in Figure 3.3. A vertex correction (see 
figure ^n| f) that contributes to this process has also been included in a second 
graph of the range 50-100 MeV. There is good agreement with the data. Our 
photodisintegration cross-section also reporoduce the correct threshold behavior 
as described by Brown Iso) . 

Another important check is ensuring that the calculation is gauge invariant. 
In electrodynamics, gauge invariance arises because Maxwell's equations remain 
unchanged under the transformation 



Approximation C: B 
Approximation D: C 
Approximation E: D 
Approximation F: E - 
Approximation G: F 
Approximation H: G 
Approximation I: H ^ 



A ^ A + VA, 



(3.65) 



where A can be any scalar function. Since the relativistic two-photon scattering 
amplitudes must have the form e^T^^e^, it follows that the condition imposed 
by gauge invariance is that 



(3.66) 



An important consequence of this is what is usually referred to as the low- 
energy theorem. Friar showed that, as the photon energy goes to zero, the 
full Compton amplitude must go as 



M 



mt 



(3.67) 



"a 
Pi 

Co 

t) 

. S 

Table 3.1: Comparison of photodisintegration cross-sections calculated in several approximations, ap is the cross-section ^ 
published by Partovi [47], while aK is from this calculation. O 
Approximation 20 MeV 80 MeV 140 MeV ^ 

O 

O 

o 
o 

Co 
CO 





ap {^ib/sr) 


aK ifJ-b/sr) 


ap [iib/sr) 


aK ifJ-b/sr) 


ap {^ib/sr) 


aK ipib/sr) 


A 


579.1 


583.3 


77.15 


80.54 


34.04 


34.56 


B 


589.2 


593.2 


83.50 


85.83 


39.77 


38.59 


C 


591.5 


595.0 


84.55 


86.86 


40.38 


38.99 


D 


591.7 


595.1 


84.85 


87.06 


40.67 


39.16 


E 


591.6 


594.7 


84.73 


86.61 


40.55 


38.70 


F 


592.1 


595.3 


85.43 


87.68 


41.31 


40.26 


G 


591.9 


594.6 


90.45 


90.22 


47.07 


43.33 


H 


588.2 


591.2 


87.52 


86.36 


44.64 


39.44 


I 


588.2 


591.2 


87.40 


86.40 


44.53 


39.52 



□ 
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Figure 3.3: Total photodisintegration cross-sections ctt vs. energy compared 
with experimental data. The range 50-100 MeV is magnified and a vertex 
correction (VC) is added. The experimental data come from |48 and p9[ . 
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This was derived assuming only that the calculation be gauge invariant. It 
is non-trivial because it involves the target (deuteron) mass; at low-energy the 
seagull amplitude depends on the proton mass : 



M 



sg 



(3.68) 



Therefore, all other terms, including meson-exchange terms, must exactly cancel 
half of the seagull amplitude at threshold. 

If we ignore meson-exchange terms, this is not as difficult as it sounds. The 
only term which survives in the low-energy limit, from equation (3.36), is 



1 



P 



(3.69) 



where we have only included the kinetic energy part of the potential energy 
term involves meson-exchange currents and will be discussed shortly. 

This can easily be shown to satisfy Friar's low-energy theorem. Since A e 
at threshold (the dipole approximation), the definition of (equations 3.24 
and 3.251 ) implies that 



e r ^ 
--•£/. 



h2 
e f 
h2 



Evaluating the commutators of equation (3.69) then yields 

^2 



M 



a4 



(3.70) 
(3.71) 

(3.72) 



which cancels half of the seagull term as required. 

Figure 3^ further emphasizes this point by demonstrating that our numerical 
calculation reproduces this result. Since the cross-section goes like the square of 
the amplitude, we expect the total cross-section at low energy to be one-fourth 
of what the cross-section would be if only the seagull term were included. The 
differential cross-section in this figure includes all terms except meson-exchange 
terms, and is seen to approach the correct limiting value. This plot shows only 
the forward cross-section, but the low energy theorem is satisfied at all angles. 

We now examine the effects of the meson-exchange currents that were pre- 
viously neglected. Since the potential energy portion of the double commutator 
term (equation 3.69) has not yet been accounted for, we suspect that this must 
cancel the explicit meson-exchange terms at low energy. This is indeed the 
case. Figure 3.5 shows the four pion-exchange diagrams which can contribute, 
according to Arenhovel ||5^ . He showed that these terms analytically cancel the 
yOPEP ^gj.^ threshold. 

These four terms and the yOP^p ^^gj-j^j calculated in Appendix I. Figure ^.6| 
shown a graph of a differential cross-section at 0° which includes only these 5 
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Figure 3.4: Forward differential cross-section vs. energy compared with the 
low-energy limit. 
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Figure 3.5: Gauge-invariant set of pion-exchange diagrams 
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Figure 3.6: Forward differential cross-section including only pion-exchange and 
yOPEP contributions vs. energy 
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terms. It is exactly zero at threshold. This is true at all angles, not just in 
the forward direction. The cross-section can be seen to have an dependence, 
which means that the amplitude goes like w^. 

This is a very non-trivial check. Gauge invariancc has been satisfied without 
explicitly demanding it in the formulation of the problem. This check, along 
with the photodisintegration comparisons, inspires confidence in the correctness 
of our calculation. 



Chapter 4 

Results and Discussion 



4.1 Effects of Various Terms in Cross-Section 

We are now ready to calculate the deuteron Compton scattering cross-sections. 
We concentrate on three energies: 49 MeV, 69 MeV (the energies of an Illinois 
experiment Q ) and 95 MeV (the energy of a new Saskatoon experiment Q ) . 

First, we study the contributions to the cross-section from the major inter- 
actions. They will be labeled as follows: 

SG: The Seagull Term (Appendix A). 

EM: All interactions arising from the multipole decomposition of the vector 
potential as described in Appendices F and G, except for recoil corrections. 
These include all electric and magnetic dipoles, as well as the double commutator 
term needed to satisfy the low-energy theorem. 



tt: All pion interactions pictured in Fig 3.5 which are needed to satisfy the 
low-energy theorem. These are calculated in Appendix J. 

a + P: The polarizability terms, calculated in Appendices B and C. Unless 
otherwise specified, the values for the polarizabilities are Op — 10.9, q;„ = 12.0, 
f3p = 3.3, (3n = 2.0, as determined in |26|| . 



These terms are displayed on the top graph in each set of graphs of Figure 4.1 
The data from |Q is displayed on all graphs at 49 and 69 MeV; no data has 
yet been pubhshed for 95 MeV. 

The seagull term alone provides a reasonable description of the data at 49 
MeV, and does no worse than all the terms together at 69 MeV. At low energy, 
the seagull term has a (1 -|- cos^ 6*) dependence. As the energy becomes higher, 
this dependence is maintained at the forward angles, while the cross-section 
gets increasingly smaller at the backward angles. It is the dominant term at all 
energies under consideration. 

The EM terms raise the seagull cross-section and become increasingly im- 
portant as the photon energy rises. The tt terms, on the other hand, do not have 
a very large effect, even at 95 MeV. These terms exactly cancel at threshold and 
continue to nearly cancel at higher energies. The amplitude of these terms has 
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an dependence, as was mentioned in section 3.6. This is borne out in the 
figure as the contribution at 95 MeV is approximately four times that at 49 
MeV. 

The polarizabihty ampHtudes also have an ui'^ dependence. They are nearly 
as large as the EM terms at energies as low as 49 MeV, and have a major effect 
on the cross-section at 95 MeV. This large effect is necessary if there is to be 
any chance of measuring the polarizabilities at these energies. 

These terms are collectively called the Major Terms (MT). They describe 
the data reasonably well without intro duci ng the smaller corrections detailed in 



the second graph of each set of Figure |4.1| . These corrections are: 



RC: The Relativistic Correction, as calculated in Appendix I. This includes 
only the first-order correction. 

CMC: Recoil Corrections (from the Center-of-Mass operator), as calculated 
in Appendix H. This also includes only the largest correction. 

VC: 7iV Vertex Corrections, such as the one pictured in Figure |T](f). These 
are calculated at the end of Appendix J. 

The RC term has an dependence. The scattering amplitude itself goes 
like u, but does not interfere with the seagull term because of its spin operator. 
It is a small effect at 49 MeV in that in cannot be detected because of the size 
of the experimental error bars. However, at 95 MeV it is important. This term 
decreases the forward cross-section and raises it at backward angles. Somewhere 
around 90 — 110° it has no effect at all. 

The CMC terms are very small at all energies here. Their presence cannot 
be detected within the error bars and do not affect the calculation. 

The VC terms tend to counteract the RC terms. Their effect is largest 
at backward angles. However, they do not increase as rapidly in energy and 
therefore become less important relative to the RC and polarizabihty terms. 
They also cancel at about the same angle as the RC terms. 

Next, we examine the EM terms in m ore detail. A plot of the cross-section 
versus energy at 9 = 0° is shown in Figure O. The following terms are included: 

SG: The Seagull Term, as above. 

DC: The Double Commutator term of equation ( 3.36 ), calculated in Ap- 
pendix F. 

EI: Electric Interactions arising from the J'^^-' operator in equation ([Gj 
calcul ated in Appendix G. This also includes the second-order terms of equa- 
tion ( [3.33| ) and ( 3.34 ). The L = 1 term is the first term in the multipole 
expansion. 

MI: Magnetic Interactions arising from the J^*^-* operator in equation ( G. 5 ), also 
calculated in Appendix G. 

The seagull term is independent of energy at the forward angle. It depends 
on q and not ujf or uji individually. The DC term is seen here to give the correct 
low-energy limit, and slowly decreases with energy. This term contains pieces of 
both magnetic and electric interactions, because current conservation was used 
to derive it. This is the reason we are not using the labels El and Ml; these 
terms in the usual sense cannot be extracted from the DC term. 
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The L = 1 terms for the EI and MI terms are added next. These terms 
are zero at low energy. The EI term for L = 1 is important at low energies, 
increasing rapidly until above 40 MeV, after which it remains approximately 
constant. The MI term with L = 1 becomes noticeable only aroimd 20 McV, 
and increases steadily. Adding the L = 2 terms produces no noticeable effect; 
we conclude that multipoles higher than L = 1 are unimportant below 100 MeV. 
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Figure 4.1: Contributions of different terms to Compton cross-section at 3 
different energies. The upper graphs include the seagull (SG), electromagnetic 
multipole interactions (Ml), pion (vr) , and polarizability terms. These terms 
are collectively called the "major terms" (MT). The second graph in each set 
includes relativistic (RC), recoil (CMC), and vertex (VC) corrections. 
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Figure 4.2: Contributions of different multipoles to the electromagnetic terms 
(MI) vs. energy. Included are the seagull (SG), double commutator (DC), 
electric multipole (EI), and magnetic multipole (MI) terms. 

4.2 Determining Polarizabilities from Deuteron 
Cross-Sections 

Having seen the size of the major contributions, we now investigate the depen- 



dence of the cross-section on the size of the neutron polarizabilities. Figures 4.3 



and 4.4 show the effects of changing the magnitudes of an and /3„ at the three 
energies we have been looking at. We first vary them one at a time, and then 
together. The term has the same (1 -|-cos^ 9) dependence as the seagull term 
at low energies, while the /3„ term goes like cos 9. This means that the magnetic 
polarizability term gives no contribution at 90° . Note that the largest effects for 
both terms are at extreme forward and backward angles, while the experimental 
data points are in the middle. 
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Figure 4.3 shows the changes to the cross-section over the range of a„ = 
12.0 ± 4.0 and /?„ — 2.0 ± 4.0. These errors are shghtly larger than the errors 
quoted in |^ but are used to illustrate the sensitivity to the size of the error 
bars. The graph of 49 MeV shows that the experimental error bars are too 
large to provide significant distinction between the values of Q!„ — 8.0 — 12.0 
and (3n = 2.0 - 6.0. 

There are only two data points at 69 MeV, and one of these is not consistent 
with the calculated cross-section. The data point at the greatest angle is also 
difficult to fit at 49 MeV. Other calculations show similar inconsistencies with 
these two points |]37|, |40| , Looking at the graphs where the polarizabilities 



both vary (Figure |4^), we see that having a lower «„ and a higher /3„ than 
the accepted values provide the best fit to these points. At 69 MeV there is 
not much else to say as the other point provides little information except that 
the calculation fits the data using an acceptable range for the values of the 
polarizabilities. 

The situation looks more promising at 95 MeV. The effects of the polar- 
izabilities are larger, and the range of ±4.0 probably will extend beyond the 
experimental error bars at angles such as 50° and 140°. 

We must keep in mind that the proton polarizabilities also have uncertainties, 
and they contribute to the deuteron cross-section just as much as the neutron. 



Another set of graphs as described above is given for ap and /3p (Figures [4.5 



and 4.6), but using only the ranges as the actual experimental error given in 
[ pT| . The same trends stand out here, most notably that a smaller value of a 
and a larger value of /3 would provide a better description of the data points at 
the back angles. 

An argument can be made that (3 has an energy dependence which has been 
thus far neglected, causing its value to be underestimated. We have already 
discussed in Section 2.2 that the A resonance gives a large contribution to the 
magnetic polarizability. For Compton scattering, this part of (3 can be written 
in terms of matrix elements as 



a" 

niA — rriN — Lu 300 — uj 



^.^,\AMMJh)l^2\i^lM^. (4.1) 



The photon energy was omitted in the static case of Section 2.1. In Figure 4.7, 
we compare the total cross-section with the static value of /3„ = 2.0 to the one 
where /3„ is a function of lo. The energy-dependent polarizability is in fact more 
consistent with the data, although the problem at the back angles still persists. 
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Figure 4.3: Effect of varying q;„ and separately in the Compton cross- 
sections at 49, 69, and 95 MeV. /?„ = 2.0 in the top graph and a„ = 12.0 in the 
bottom. 
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Figure 4.4: Effect of varying both a„ and /3„ in the Compton cross-sections at 
49, 69, and 95 MeV. 
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Figure 4.5: Effect of varying Up and Pp separately in the Compton cross- 
sections at 49, 69, and 95 MeV. /3„ = 3.3 in the top graph and a„ = 10.9 in the 
bottom. 
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Figure 4.6: Effect of varying both a.^, and /3p in the Compton cross-sections at 
49, 69, and 95 MeV. 
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Figure 4.7: Comparison of energy-dependent vs. static P in the Compton cross- 
sections at 49, 69, and 95 MeV. 
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4.3 Sources of Error 

We would like to study the effects of some of the smaller terms more closely. The 
results presented here basically reinforce the general discussion of Section 4.1. 
For each term, we show a graph comparing the full cross-section with a cross- 
section that has that effect subtracted out. We can then estimate the absolute 
error in the value q;„ (we fix /3„ = 2.0 in this section) caused by neglecting 



various effects. These are listed in Table 4.1 and shown in the figures 



First, we examine the gauge invariant set of pion terms. Examining Ta- 



ble [4.1| , we see that the value a„ would be slightly underestimated at the forward 
angles and overestimated at backward angles, and would not be affected around 
120 — 140°. The magnitude of this error is about the same as the experimental 
error in the polarizability, so these effects would certainly need to be taken into 
account to make an accurate determination of a„ in this type of experiment. 
The effect is same at all energies since both the tt terms and polarizability terms 



have an dependence. This is shown in Figure 4.8. 

The recoil (CMC) corrections are difficult to see on the graphs (Figure [ 
and this is also shown in the table. The largest errors are actually at the lowest 
energies, indicating that these terms are of order lu^. The small values for 
Aa„ at 95 MeV (especially at forward angles) indicate that these terms can be 
neglected here. 



The relativistic corrections (RC) (Figure 4.10) are not a "small" correction: 
even at 49 MeV they would have to be included in a calculation to make an 
accurate determination of a„. The amount of error that would be introduced 
by neglecting these terms is comparable to a„ itself at some angles. They also 
scale like so the energy dependence should not be very strong. This is the 
case at backward angles, but at forward angles this effect increases with energy, 
perhaps due to interference with another term. Again, this effect is smallest 
around 110° at each energy. 

We turn to the vertex corrections. These terms would also need to be in- 
cluded for an accurate estimate of the polarizability. Their effect is actually 
greater than the effects of all of the other pion terms combined. The errors de- 
crease as energy increases, but even at 95 MeV the errors, particularly at middle 
angles, are at least as large as previous errors in polarizability determinations. 
At backward angles the effect seems to be independent of energy. The graphs 



are shown in Figure 4.11 



Finally, we make an estimate of the effect of the uncertainty in the ttN 
coupling constant /. We have been using the value p /[Tic) = 0.075, but the 
uncertainty in this makes a value of 0.08 also reasonable. The effect of this is 
not large, as seen in Figure 4.12| ), but it can be seen to have a minor effect at 



higher energies in the middle angles. 
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Table 4.1: Absolute error in electric polarizability generated by omitting var- 
ious terms - pion (tt), recoil (CMC), relativistic (RC), and vertex (VC) correc- 
tions - at 49, 69, and 95 McV as a function of angle. The dependence of the 
cross-section on the value of is also examined. 
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Figure 4.8: Effect of neglecting pion terms on the full differential cross-section 
(FC) at 49, 69, and 95 MeV. 
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Figure 4.9: Effect of neglecting recoil terms (CMC) on the full differential 
cross-section (FC) at 49, 69, and 95 MeV. 
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Figure 4.10: Effect of neglecting relativistic terms (RC) on the full differential 
cross-section (FC) at 49, 69, and 95 MeV. 
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Figure 4.11: Effect of neglecting vertex corrections (VC) on the full differential 
cross-section (FC) at 49, 69, and 95 MeV. 
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Figure 4.12: Effect of changing on the full differential cross-section at 49, 
69, and 95 MeV. 
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4.4 Tensor-Polarized Deuteron Compton Scat- 
tering 

Another interesting observable is the tensor-polarized deuteron Compton scat- 
tering cross-section. It is defined by 









't 4 



d<7 , , da , , da , 



(4.2) 



This observable has not been well-studied; a recent calculation in effective field 
theory can be found in 1 53 . Although this cross-section does indeed have a slight 



dependence on the polarizabilities, we do not expect this to give a measurement 
of a and /3. However, it is easy to calculate within our formalism, and we hope 
that our p redic tions will be of use to others interested in this quantity. 

Figure 4.13 shows the total polarized cross-section for two different energies, 
both with and without the polarizabilities. These results are similar to the 
unpolarized case in that the effect of a and /3 is greater at forward angles and 
at higher energies. However, the absolute magnitude of this effect is at most 
10% of its magnitude in the unpolarized case. Moreover, the polarized cross- 
section is much more sensitive to individual terms, as we will show shortly. It 
does not seem that this quantity could produce an accurate measurement of the 
polarizabilities, but since no experiment has been performed yet, we must wait 
for this to be done to say anything for certain. 
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In an effort to be more consistent with the calculation of |53 , we now try 
to isolate the contributions of the magnetic interactions and the pion-exchange 
terms at one of these energies. The pion-exchange terms include any pions that 
give rise to the tensor force (and the I — 2 state) in the deuteron. However, we 
will include only one pion exchange in the pion terms and no pion exchanges in 
the magnetic terms. Therefore, the D-st&te deuteron wavefunction U2{r) will be 
set to zero in all terms except for the Thomson term; the D-state contributions 
to this term will be grouped with the pion-exchange terms. 

The results for the magnetic interactions (MI) are shown in Figure 4.14 . 
These interactions include all of the terms that are calculated in Appendices 
F and G, as well as the seagull term (as defined in Appendix A) without the 
deuteron D-state. In addition, we show the effect of including the deuteron 
i)-state in the calculation of the magnetic terms (but not the seagull term). 
This is a non-negligible effect. Finally, the relativistic correction (RC) is added 
(including the D-state). Since this is a spin-orbit effect, we group this with 
the magnetic interactions. Both the relativistic term and the deuteron D-states 
have a large effect in the calculation; this observable is very sensitive to these 
individual effec ts. 

Figure 4.15 shows the pion terms. This includes all terms in Appendix J, 
both the gauge invariant set and the vertex corrections, along with the D-state 
contributions of the seagull term. We also show the effect of including the 
deuteron I?-state in the pion terms of Appendix J. This effect is not a large as 
for the magnetic interaction case. 

Finally, we would like to estimate the magnitude of the interference of the 



pion and magnetic terms with each other. The solid curves of Figures 4.14 



and 4.15 show the interference of the magnetic and pion terms with the S- 
state seagull term (as well as with themselves), since the S'-state seagull term 
alone gives no contribution to the tensor-polarized cross-section. The dashed 
curve in Figure 4.16| shows the sum of these two curves. In addition, the full 
calculation (not including polarizabilities) is shown as the solid curve. The 
difference between these curves is the interference of the magnetic and pion 
interactions with each other. This is seen to be a minor effect, but is largest at 
forward and backward angles. 
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Figure 4.13: Tensor-polarized deuteron Compton scattering cross-section, with 
and without polarizabihties. The values of an = 12.0, /3„ = 2.0, ap = 10.9, and 
(3p = 3.3 are used in the solid curve. Both curves include all other interactions 
described in the Appendices 
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Figure 4.14: Tensor-polarized deuteron Compton scattering cross-section, in- 
cluding only magnetic interactions(MI). The dotted curve contains all interac- 
tions in Appendices F and G, as well as the S'-state seagull term, but only the 
contributions from the deuteron 5-state are included. These contributions are 
added in the dashed curve (but still only the 5-statc is included in the seagull 
term). Relativistic corrections (RC) are added in the solid curve. 
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Figure 4.15: Tensor-polarized deuteron Compton scattering cross-section, in- 
cluding only pion intc;ractions (tt). These consist of the gauge- invariant set of 
terms described in Section 3.6, as well as the vertex corrections and contribu- 
tions to the seagull term from the deuteron Z?-state. The dashed curve includes 
only the deuteron S'-state in the gauge-invariant pion terms and the vertex 
corrections, while the solid curve includes both the S- and £>-states. 
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Figure 4.16: Effect of interference terms between the pion and magnetic terms 
in tensor-polarized deuteron Compton scattering cross-section. The dashed 
curve is the sum of the total cross-sections from the pion and magnetic terms 
individually, while the solid curve is the total cross-section including both the 
pion and magnetic terms (but not the polarizabilities). All deuteron £>-states 
are included. The difference between the curves is the interference between the 
pion and magnetic interactions. 
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Conclusion 



We have presented a calculation of deuteron Compton scattering valid for ener- 
gies less than 100 MeV. All terms that we believe to be important are included in 
this formulation. The seagull, lowest-order multipoles, and polarizability effects 
dominate the cross-section. Relativistic corrections must also be incorporated in 
the calculation. Although smaller than the preceding terms, pion-exchange and 
vertex corrections are also needed to accurately determine the polarizabilities. 
We conclude that recoil effects and multipoles of second-order and higher are 
negligible in this energy range. 

This is the most extensive study so far of the effects of various terms in the 
cross-section on the value of «„. However, certain corrections which are believed 
to be small have been omitted. No contributions from the A resonance have 
been explicitly included, although the magnetic polarizability itself is thought 
to contain a large paramagnetic contribution from the A. The A resonance has 
been shown to have a small effect on deuteron photodisintegration at energies 
less than 100 MeV [|4|. In addition, considering the rather small contributions 
from TT-exchange currents, heavier mesons such as the p have also been neglected. 
Higher-order relativistic corrections have not been included, but may need to 
be investigated further at energies near 100 MeV because of the large effects of 
the leading-order terms there. 

The only full calculation of this type that has been published to date was 
by Weyrauch Q . There are several important differences between his cal- 
culation and the present one. One of these is the NN potential used. He works 
with a separable NN T matrix, while here we have found the Green's function 
for an intermediate NN state interacting via the Reid93 potential. He also 
assumes meson-exchange currents and relativistic corrections to be negligible. 
Probably the most important difference, however, is in gauge invariance. Gauge 
invariance is broken in the separable potential model and needs to be artifi- 
cially restored. We have made no such assumptions here. Gauge invariance is 
a natural consequence of including all of the appropriate terms. 

The cross-section predicted by Weyrauch is higher than the experimental 
data. However, results published recently Arenhovel and Levchuk and L'vov 



82 



CHAPTER 5. CONCLUSION 



83 



[ p8| p9| are in agreement with the data. They have also used a diagrammatic 
approach and have included realistic potentials and meson-exchange currents, 
but the full details of their calculations have not yet been made available. Lastly, 



a new calculation using effective field theories also agrees with the data |40 



We return to the basic question of the feasibility of making an accurate 
determination of the neutron polarizability from a deuteron Compton scattering 
experiment. The available data at 49 and 69 MeV show that the ranges q;„ = 
12.0 ± 4.0 and /?„ = 2.0 ± 4.0 are reasonable, but a more accurate estimate 
is difficult since the polarizability terms are small relative to the experimental 
error bars. The uncertainties in the proton polarizabilities must also be taken 
into account. These problems are compounded by the data point at the greatest 
angle at each energy, which we find difficult to describe (other calculations also 
find disagreement with this point). The closest fit to the all of the experimental 
data comes from lowering the value of a„ and raising Throwing out the 
data point in question allows a closer fit, but there would still be a wide range 
of reasonable values for the polarizability. 

Perhaps the most promising estimates of the polarizability will come from 
data at higher energies. Assuming that the experimental error bars are not 
significantly larger, we should find that the range of polarizabilities that fit 
the data is smaller. Sensitivity to both a and /3 is largest at the forward and 
backward angles, while the cross-section at 90° is the least sensitive to the 
polarizabilities. Many of the "small" corrections, however, are particularly small 
at 120 — 140°, and the cross-section is still very sensitive to the polarizabilities at 
these angles. The proton polarizabilities, unfortunately, will also make a larger 
contribution to the cross-section, and there is no way to determine the neutron 
polarizabilities more accurately than those of the proton in this experiment. 
This energy is at the limit of the validity of the low-energy expansion, which 
may introduce errors as well. In any case, studies at higher energies should 
provide additional information on the polarizabilities. 
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Appendix A 

Seagull Term 



The interaction Hamiltonian is 

j—n,p ^ 

since the neutron has no charge. Using 



^(^"^•) = -^ E ^J^[akexe''■'^+alile-'^^■% (A.2) 

fc,A=±l 

we obtain the transition matrix 

T/^« = (d/,P/,7/|i?^^|d.,P.,7.) (A.3) 

J nip V ^i^iOJf ' 

{ df, Pf I e-i«-^^e-^«-^ \f,R){f,R\di,Pi), (A.5) 

whereg= kf-ki, r = Xp-Xn , -R = {xp + Xn)/2, p = pp-pn, and P = Pp+Pn. 
The notation 7}j means that T is a function of the final kinematic variables 
kf, Mf, A/ as well as of the corresponding initial ones. Assuming the center-of- 
mass wavefunctions for the deuteron states to be plane waves, the integral over 
R can be performed to obtain 

(A.6) 

We now define the scattering amplitude M.^^ by 



_ 1 2-^% 5{Pf + kf-Pi-ki) 
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and so 

;=0,2('=0,2"' P 

(A.8) 

The deuteron wavefunction can be written as 

(r|^.) = ^3^r(r) + ^3^r(^), (A.9) 

where Uq, U2 are the radial wavefunctions for the 5, D states, respectively, 
and y^'^{r) = { r \ LSJM ) is the combination of the angular and spin 
wavefunctions. Inserting the partial wave decomposition of the exponential, 

oo L 

e-^'"=E E M-i)''jL{^)Y£MmLM{f), (A.10) 

L=OM=-L 

and the wavefunctions into the matrix element gives 

WLM-' ^ 

u,{r)uAr){ I'llMf \ Ylm \ HIM, ). (A.ll) 
We can use the Wigner-Eckart theorem on the Ylm matrix clement to get 

{i'nMf\Yr.M\inM,) = {-iy-''f(^_l^^ ^ ^J(m||yi||zii). 

(A.12) 

This reduced matrix element can be further simplified. First, we rewrite Ylm 
as a tensor product [Yj, 1]lm, where 1 is an operator in spin space and where 
the tensor product is defined by 

ni»i2 

(A.13) 

We then separate the spin and orbital angular momentum parts using the rela- 
tionship 

( jihJW [Tk, ® Tfejfe \\n'j2'j' ) = (A.14) 

V(2i + l)(2j' + l)(2fc + l) I n' 5 f \ { JiW Tk, )( J2II \\j2' ). 

ki k2 k 



This gives 

I' 1 1 

-Mf M Mi 



L L 
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= (-1) 



3V2L + 1{ mil 110(111 1111 
1 L 1 
-Mf M Mi 

3( /'II Yl \\l ), 



l-M 



1 I' 1 
I 1 I 



where the 9j symbol has been simpUfied using 




V'(2e+1)(2/ + 1) \ d c f 



(A.15) 

}^ 

(A.16) 
(A.17) 



and ( 1|| 1 ||1 ) = \/3 has been found from the Wigner-Eckart theorem. The 

formulas for several useful reduced matrix elements, including ( /'|1 Yl \\l ), can 
be be found in Appendix K. Putting everything together into Mff', and then 
removing the formal sum over M by substituting M = Mf — Mi, we obtain 



tft^r) 



(A.18) 



'IL 



yIm.-m, {q)VM'2i' + i)(2i + 1)(2/ + 1) <j } [ I 



I L I' 




1 L 

-Mf Mf - Mi M; 



These sums can be evaluated explicitly since there are only 5 triplets {l',L,l) 
which give a non-zero contribution. Defining 



li'^ dru,{r)j,{^)u,{r), 



we can write the final result as 



M 



SG 
fi 



rrir, 



(A.19) 



(A.20) 



1 1 

-Mf Mf 

1 2 1 

-Mf Mf - Mi M, 



t22 ' 



Appendix B 

Electric Polarizability Term 



The interaction Hamiltonian is 



aj dA{xj) dA{xj) 



2 dt 



dt 



We insert equation (A. 2) for the vector potential into 7^^ to get 



an (eA.e''^'-^" 

LOfLOi 2TTh 



2( df,Pf\ 



I "p ( 



-ikf -Xj 



After integrating out the center-of-mass states, this becomes 

M% = -LVfij.iil^^ ■ exj J2 J2 ^ 

1=0,2 i'=a,2 



(B.l) 



(B.2) 



(B.4) 



dV( I'UMfrf I (ape-59-*='+ a„e59-*') | f){r\lllM,ri ). 



We can expand the exponentials into partial waves using equation (A. 10) and 
insert the deuteron wavefunctions (equation [A.9|) to get 



oo L 



I'l L=0 M=-L''^ 



uAr)j,{^)u,{r) X (B.5) 
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This is identical to equation ( A. 11 ) if we replace ^ by —ujfUJi \ap + (— l)^a„] , 
so we can im medi ately write down the final answer by making this substitution 
in equation (A.2C). Thus, 



M% - -c^/c^,Vl27r(e^^ •eAj(-l)'"*'^K + a„) X 



-Mf M 

1 2 1 

-Mf Mf-Mi M, 



t20 



(B.6) 



t22 "1 
— ) 



The factor (—1)'^ has been omitted since L must be even. 



Appendix C 

Magnetic Polarizability 
Term 



The interaction Hamiltonian is 



WjxA{xj) ■ W,xA{x,) . (C.l) 



j=n,p 



Using this to calculate TP- gives 



and so 



(C.2) 



dV( I (/3pe-*«-^~ + /3„e*«"-^^) | f )( f| ). (C.3) 



This is the same as equation (B.4), but with a replaced by (3, and with a differ ent 
dot product. Therefore, we can make these substitutions in equation (B.6) to 
get 



-ujfUJi\/l2i: 



(fc. X ex) ■ [kf X elj] (-l)i-A^/(/3p + /?„) 
Yoomio° + Io')5M,.M.+ 



1 1 

-Mf Mf 

1 2 1 

-Mf Mf - M, M, 



t22 

02 , t20 ^2 



V2* 
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Appendix D 

Phase Space Factors 



To calculate the differential cross section, we use Fermi's golden rule 

final 7 final d 

1/V (27r)3 (27r) 

initial 7 flux 



27r?i 1 



(D.2) 



;i l/V (27r)3 (27r)3 ^/tUyij- 1/2 
= -^d{AE)d^kfd^pfS{Ap)\Mff , (D.3) 

where we have used the definition of Mf^ (equation |A.7] ), as well as |(5(AP)|2 = 

T/(27r)3^(AP). Performing the integration over pf removes one of the (5-functions; 
da becomes 

dafi = n ^^^ ~ ^K /dujfdnf\Mff , (DA) 

OJfUJi 

and so 

where E{E') is the total energy of the initial(final) deuteron and photon com- 
bined. The partial derivative is calculated in the lab frame. We define 4-vectors 
pi = (TOd,0), Pf = {Ef,pf), ki = {u>i,ki), and fc/ = (tj/,fc/). By energy- 
momentum conservation, pi + ki = pf + kf. Moving kf to the left-hand side 
and then squaring both sides gives 

h-Pz-kf -pi-kf - ki = 0. (D.6) 
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This simplifies to 



rrid + UJi — uJi cos ( 



(D.7) 



where 9 is the angle between fcj and ki. Next, we use energy conservation 
(Ef — LUi + rrid — LOf) to write 

Ef = \J cjj2 + -I- 71^2 _ 2uofU)i COS0, (D.8) 

where we have again used equation (pD.Tj ) . Now we take the derivative: 

(D.9) 

1 



dE' 
dujf 



d{Ef+ojf) 

dujf 
ujf — UJi cos t 



Ef 



This has been simplified using equations (D.7) and (D.< 
into equation (D.5), we see that 



da_\ 

dn) 



Ef 



M 



fi I ; 



(D.IO) 
(D.ll) 

Plugging this back 
(D.12) 



where M J:^ — Mff' + MJ^ + M'^^ + ■ ■ ■ . If we are interested in computing the 
spin-averaged cross section, we sum over the final and average over the initial 
polarizations to get 



all 6 -^-^ V LUi J rriri ■' 



MfMiXfXi 



(D.13) 



Appendix E 

Multipole Expansion 
Derivation 



It is convenient to obtain the multipole expansion for i\e'^'^'^ , where e\ is a spher- 
ical vector. This derivation is similar to the one found in some additional 
details are included here. 

We set k to be in the direction of the z-axis for now, which means that 
ex becomes f\, the unit vector in the spherical basis, with A = ±1. Defining 
a;= Ifcl, 



ik-r 



k—z 



^^Y.T. ^^^'M^r)Y;^ik)YUr) (E.l) 
fx v/47r(2/ + {u:r)Ym (r). 



(E.2) 



where we have used equation (A. 10) and Yi„i{z) = (5mOY ^37^- "^^^ define 
the vector spherical harmonic Tj l M{f) by 



-JLM 



L J 

-V M + V ~M 



By rearranging this to get 



f-vYLM+vif). 

(E.3) 



Y^(f)f-x = ^(-1)-'+^^+i^/2JTT _ 



1 I J 

-A m -M 



TjzA/(f), (E.4) 



we can simplify equation (E.2): 

Ak-r\ 



I k—z 



(E.5) 



j=\i-i\ 1=0 



E E(-l)'^'"'*V47r(2/ + l)(2J+l)f I J ^^\Mujr)fjix{f) 
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Explicitly writing out the J sum yields 

Ak-r\ 



1=0 ^ 



V21 + 1 

oo 

= i^\f^ji {lot) 



1 / / 

A -A 



1 I /-I 
A -A 

1 / l + l 
A -A 



1=0 



I - 1 



-Ti-iix{fy 



(E.6) 

Ti+ii\{f) 

(E.7) 



after evaluating the 3j symbols. Next, we manipulate the summation index in 
order to rearrange this into a form where all of the vector spherical harmonics 
have the same total angular momentum. We call this new index L, and let 
L = / — 1 in the first term of equation (E.7), i = / in the second term, and 
L = / + 1 in the third term. Thus, 



Cab 



ik-r 



k—z 



L = l 



L 



2L + 1 



iL+i{ujr)TLL+i\{r)- (E.; 



\j L{ujr)TL L \{f ) - i 



L + l 



2L + 1 



Now we would like to return k to its original orientation, where it has spherical 
coordinates {Lo,'d,(f). This can be attained by rotating r by the Euler angles 
(0, — t?, —if) by means of the Wigner d-functions L)^^|(0, — — <p), so that 

oc L I 



CAE 



L = l M=-L 



L 



-j,^,{^r)TLL+iM[f)- (E.9) 



\jL{(^r)TLLM{f) - i 



L + l 
2L+1 



2L + V 

J.-iK)fiL_iM(r) \ D^iil{0,-d,-cp). 



Before simplifying this further, we need to show that 



1 



VL(L + T) 



1 



ujfj^{ujr)YLM{r) + -V,. 1 + r— ji,(wr)lLM(r) 



dr 



L + l 



2L + 1 



3L-i{^r)TLL-iM{-f) - 



2L + 1 



jL+i{ujr)fLL+iM{f){E.W) 



and 



LYLM{r)jL{^r) = 3L{i^r)TL l M{r). 



(E.ll) 
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We start with the first term on the left-hand side (LHSl) of equation (E.IO) 
LHSl = 



^ J j,K)W^(-l)'-yi.(f)r,.r^M(r)(E.12) 



where we have used the decomposition of an arbitrary vector V into spherical 
components, 

V= (-l)"Kr-., (E.13) 



along with 



yyr.(n 



Next, the spherical harmonics are combined using 



Y Y-, , 



j+j' J 
J=\j-j'\M=-J 



M 



(2j + l)(2/ + l)(2J+l) 
An 



J f J\f j f J 
/ I TO to' -M 



YjM 



(E.14) 

(E.15) 
(E.16) 



to get 
LHSl 



l+L J 



E E E 



j,(fcr)(-l)-+^V(2L + l)(2J+l) 



I L J\ f I L J 







V M -M 



(E.17) 



Looking back at equation (E.4), we see that the spherical harmonic can be 
replaced by a vector spherical harmonic: 



LHSl 



l+L 1+J 

E E E (-i)^-W(^ 

u=-l,0,l j=|i_L| .7' = |1-J| 



(2L + 1)(2J+1)(2J' + 1) 



1 L J 




T 



\j'JM'if) E ( ^ 



MM' 



1 L J 

M -M 



L{L + 1) 



This can be simplified because of the orthogonality of the 3j symbols. 



E 



Jl J2 ]3 
mi 7712 TO3 / \ nil 7712 TO3 



jl j2 j's \ _ Sj^j^'S; 



i'Omsms' 



2h + 1 



(E.19) 



to yield 

LHSl = 



2J+1 fl L J 



L{L + 1)\ J ^-^^^ 



f^,;^,(f)(.E.20) 
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J can only be equal to L + 1 or L — 1 . Inserting this into equation ( E.20| ) and 
then writing out the 3j symbols algebraically, we obtain 



LHSl = (-l)^wji(wr) 



' 2L + 3 / I L L + 1 
L{L + 1)\0 



T, 



L L+1 M 



2L- 1 



1 i L- 1 
+ V 



TLL-lM{f) 



TLL+lAli f) ^ 



(E.21) 



(E.22) 



The second half of the left side of equation (E.IO) can be simplified with repeated 
applications of the gradient formula (equation K.13|) and the recursion relations 
for spherical Bessel functions (equations K.ll and pK!.12) to get 



LHS2 : 



Tl 



L+IM 



(f) 



V^(2i + 1) 
ujr- — ] jL{ujr) 



2L+ 1 



2L + 1 



TLL-lM{f) 

v/(i + l)(2L + l) 



L{L + r 



2L + 1 



or ^ J^+iK + ^ - j^iior) 



Combining this with equation (E.22) and then using the Bessel function recur- 
sion relations again gives 



LHSl + LHS2 = 

ThL+lMif) 



(E.24) 



^L{2L + l) 



l2 

ujr 



L2 L{L + 1) . 



2L+1 
v/(i + l)(2L + l) 



2L + 1 

L{L + 1) . 
2L + 1 ^ 



L 



2L + 1 



jL + i[l^r)TLL+lM{f) + 



L + 1 
2L+ 1 



3^_,{Lor)TLL-iMm-2b) 



This is identical to equation (E.IO). To prove equation ( p. 11 ), we first write o ut 
the left-hand side in terms of its spherical components using equation ( E.13| ): 



.LYLMif) = J2 



i-ir 



^/L{L + 1 
Using the fact that |||, p. 122] 

(_l)-Af-<^+i+i^^(2. + 1)(2L + 1 



L,f_,YLM{f). (E.26) 



LM 



L 1 

M V 



L 



(E.27) 
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this becomes 



l=LY,M{f) = ^(-l)i+^-*^V2rTT 
1^ + ^) 



L 1 L 

M V -M - V 



With the help of equation (E.3) we can immediately write 

1 



^L{L + l) 



(E.28) 
(E.29) 



which is what we w ere trying to esta blish. After substituting equation ( E.29 ) 
and equation ( E.25 ) in equation ( E^ ) , we get the final result for the multipole 
expansion: 



,L 27r(2L+l) 



L=l M=-L 



L{L + l) 



(E.30) 



d 



V.r ( 1 + r— ] iL{^r)YLM{r) - iujfj,:,{ujr)YLM{f) - \LYLM{r)ji^{Ljr) 



Appendix F 

Largest Dispersive Terms 



The remaining terms will be calculated using second-order perturbation theory: 
df,Pf, 7/ I H^'^' \C,Pc){C,Pc\ H'^' I d„ P.,7») , ^) 




hwr + Ed, - Ec - Pl/2md + ie 

{df,Pf,-ff\H"^'\C,Pc,j^,jf){C,Pc,l^,jf\H^^'\d.,P,,^,) \ 
-hu;-f + Ed,-'Ec-Pc/2md + ie J" 

Here C denotes the internal quantum numbers of the intermediate np state, and 
is shorthand for all sums and integrals which are needed to describe this 
complete set of states except for the center-of-mass states. These are written 
separately as Pc- At this point, it is most convenient to use i?™* = — / J(^) • 
and then expand A into multipoles according to equation ( E.3C| ). There 



are 3 terms in this expansion; therefore, equation ( |F.l| ) contains a total of 18 
terms. The largest contributions at low energies (which we will call T^") arise 



from the gradient operator in equation (E.3C). These are the leading-order 
contributions; the other terms are of order Lo/rriN smaller. This gradient term 
includes the case where there is an El interaction at both "fN vertices. 
With this in mind, we define ^i{f) and ^f{r) by 



1 Ph^^ A fi) , /27r(2i + l) , , 



1+?-^ ]JL{uJ^r)YLM{f) 




' 1 + 7-— ]L'{uJsr)YL'M-[r), 



L'{L' + l) \ dr 



100 



APPENDIX F. LARGEST DISPERSIVE TERMS 



101 



which means that 



1 2TTh. 



iki -r 



l2Trh 



V.$/(r) 



(F.4) 
(F.5) 



We have chosen a coordinate system where ki hes on the z-axis and where fc/ 
has spherical coordinates (c j/-, 9 , (p). Using equations (F\4) and (F\5) in the 
uncrossed term of equation (F.l) gives 

vr" -J2(df,Pf\J /(f) • ^e'^fiiv'e \c,Pc)x (F.6) 



C,Pc 

1 

huoi + Ed, ~Ec- Pll^md + is * 
Integrating by parts and then using current conservation, 



j=n,p 



gives 



'c)x 



(F.7) 
(F.8) 

(F.9) 



C,Pc 



1 



n'ihuj^ + Ed, -Ec- P^/2md + is) 
Performing the integrations yields 



C,Pc\ [H,piO]MOd'ad^,P^). 



^{df,Pf\[H,ie/h)^fixp)]\C,Pc 

C.Pc 

1 



(F.IO) 

C,Pc I [H,{e/n)MS!p)]\d^,P^). 



huJi + Ed, -Ec- Pc/2md + ie 
It is easiest to evaluate 7^^ by writing the commutators in terms of f and R 



[H,{e/n)Mxp)] = 



2mp 2m„ 



+ V,{e/h)Mxp) 



-V,.$,(r/2)e''='--^ P 



^ + y,(e/n)$,(f/2) 

Too 



Jki-R 



(F.ll) 



(F.12) 
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and 



[H, {e/h)^fixp)] 



le 



-ik f R 



(F.13) 



where equations (F.4) and (F.5) have been used, along with p = (pp—pn) /2 and 
P = Pp+Pn- The terms containing P are recoil corrections and will be calculated 
in Appendix H. Defining the dimensionless functions $i = {e/h)^i{r/2) and 
= {e/h)^f{r/2), along with H"p = ^+ V,we get 



E 

C,Pc 



df,Pf\ [H"P,^f]e-'''f^ \C,Pc){C,Pc\ [i?"P,$,]e*-« | d„P, 



huj, + Ed, -Ec- P^/2md + i£ 

(F.14) 

After evaluating the commutators and integrating out the center-of-mass mo- 
tion, this becomes 



<5(P, + fc. - Pf -kf)\]^{df\ [i?"^ I 



df\^f[H-rn\d^ 



TlLJi 



2md 



df I I d. 



(huji) 



2l2 



2m„ 



E 



I I C )(C|<1>. Id.) 



(F.15) 



2md 



Similar treatment on the crossed term of equation (B\l) yields 

1 



i(d/|<lai7"f,$/] Id, 



hujf + 



hujf 



2md 



E' 



2md 

d/|$, |C)(C|$/|d. 



2md 



d/ I I d, ) + 

(F.16) 



■ Edi - Ec + is 



We now add equations (F.15) and ( F.16[ ) to get 

r/, = 5{P, + h- Pf - kf)[Tfl + Tff + Tff + TfX 

where 

(TicOnV (d/l«'/|C)(C|$,|d,) 



hujf 



2md 
2md 



E 

c 

2 

E 



2nid 



-Ec + ie 



{dj\^,\C){C\^f\d,) 



C ~^f--2it 



Ed, - Ec + ie 



(F.17) 

(F.18) 
(F.19) 
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/7-a4 



f H — nuji 



2md 



df I I d, ), 



(F.20) 
(F.21) 



T^'*! will be calculated first. We write out $i and $ / expl icitly (see equations 
and F\3), and then convert T^^ to A^fj using equation ( A. 7 ) to get 



M 



2md 



C LM L'M' 



, / 27r(2L' + l) 



(-1) 



1 



hujj 



i^-^y ^ J. -IT-— I ^ ^ I (0,0,0)^ X 



/27r(2i + 1) 



^ + r^]jL{^)YLM{f) \d, ). 



(F.22) 



The sum over C can be explicitly divided into radial, angular, and spin sums; 
the radial sum can be collapsed. Removing the complete sets of angular and 
spin states (which reduces d'^r to dr) yields 

2 



M 



27re^ 



{huji){hujf) 



2md 



EE E E E (-1)^' 



-Xf-L+L\ 



1=0,2 l'=0,2 LcScJc^Ic LM L'M' 



r^4r r",r'B[-;i (0, 0, OD^l,, (0, 0^^^L^ 



2ma 



1 



[l + I'llMf I Yl'm' I LcScJcMc ) 



LcScJcMc I Ylm I ^IIM, 



1) 

I r ) X 
(F.23) 



{huji 



Defining Eq = 
can rearrange to get 

Mf ^ 



Eh, where Eh is the deuteron binding energy, we 

{hi,,? ^ 



2md 



E E EE(-i)^' 

W LcScJcMc LM L'M' 



-Xf-L+L'^ 



i?i^\^(O,O,O)i?i^,l^/O,-0,-^ 



' (2L+ 1)(2L^ + 1) 
LL'{L + l){L' + 1) 



dr dr' r r'u,{r)u,{r')i,,,(^)G,^{r,r'-Eo)^j,,(^) x 
I'llMf I Yl'M' I LcScJcMc )( LcScJcMc \ Ylm \ HIM, >(F.24) 
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where G^^ (r, r'; £^0) 

by 



r ). We have also defined a function ipL 



'^l{x) = jL{x) +x—jL{x). 

ax 



(F.25) 



The M sums are formal and are removed to produce the final result 

T 2 00 



{hLOi){flU) f) 



2md 



flOJj 



00 l+L Jc+Sc 

H H Ex 

/=0,2/'=0,2 L.L' = l Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ 



(2L + 1)(2L' + 1) / 1 



L' 



Jc 



Y LL'iL + 1) (L' + 1) V -M f Mf-M,- \ X, + M, 
( I ) ( Yl' WLcScJc )( LcScJcW Yl \\in ) x 

(F.26) 



dr dr' r r'u,{r)u,{r')4^,.(^)G,^{ry ■Eo)i,,(^). 



Equation ( F.19 ) can be handled similarly. The result is 

)n 2 00 



M'f: = 



27re2 L {hcjff 



.DO 1+L Jc+Sc 

E E E E E E 

l=0,2l'=0,2 L,L'=1 Sc=0,l .Jc = \l-L\ Lc = \Jc-Sc\ 



/ (2L + 1)(2£' + 1) 
LL'{L + 1){U + 1) 



1 L 



Jc 



where E'q = 



Jc L' 1 

-Mf Xi ~Xi + Mf )\Xi-Mf Mf-Mi-Xi Mi 

I'llW Yl WLcScJc )( LcScJcW Yl' ) x 
dr dr' r r'u,{r)u,{r')i^,,i^)G,Jr,r';E'„)M'^), 



(F.27) 



Equation ( |F.20|) contains no intermediate states and therefore no Green's 
function is needed. Using the definitions of and $/ from equations (F.2) and 
we get 



2^e2 



{hujf){huji 



2m<j 



2nid 



EEE(-i)""'^^"^''x 



IV LML'M' 



(2L + 1)(2L' + 1) 



I'llMf I YlmYl'm' I lllM, ) X 
dru,(r)M,/(r)V'j,'(^)-!/;i.(^Y^). 



(F.28) 
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The spherical harmonics are combined according to equation ( K.10| ) to produce: 



Irud Zrud 

(_i)L'-A,+M^(L)^^ (0, 0, mfplx, (0, -e. - 



IV LML'M'Im 

I 2L + 1 
LU{L + l){L' + l) 





L' 




I 








{ I'llMf I Y^j^, I lllM. ) / dr u,[r)uAr)^.'(^)^.(^). (F.29) 



Finally, the M sums are removed to yield the answer 

a3 e2 0F 



M 



{hiUf){huj,) ; 



/jw y H rtCJ- 



2md 



2md 



L+L oo oc 

E E E EEx 

;=0,2 ;'=0,2 L=|L-L'| i = l i' = l 



2L + 1 





L' 




V 








L 
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The Hamiltonians in the double commutators in the T^^ term (equation F.21 ) 
will be split into kinetic and potential energy parts. Inserting a complete set of 
r states into the kinetic energy term, we see that 
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{r\d, )M?A, (0, 0, 0)D'^,>_^^ (0, -0, -</>). 
Evaluating the double commutators yields 



(F.31) 
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Next, using the identity Vf -Vg^ ^[V^C/y) - fV^g - gV^f], we get 
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(F.33) 
(F.34) 



each of the three terms inside the bracket in equation (F.33) can be simphfied. 
Thus, 
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where equation ( K.IO ) has been used to combine the spherical harmonics. After 
some algebra this becomes 
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r- ■ - 2 2 

We now remove the formal sums over M to obtain the final result: 
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The potential energy part of the double commutator term, T^^'^^, will be cal- 
culated in Appendix J using the one pion exchange potential. 



Appendix G 

Other Dispersive Terms 



The next set of terms (that we call T^^) to be investigated are those in equa- 
tion (F\l) where the substitution A = V<I> is made only once. This includes 
the situation where there is an El interaction at one vertex, and a different 
type of interaction at the other. There are four types of terms: 
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(G.l) 



Using the same steps that were used to derive equation ( F.15 ) (again ignoring 
the P terms in equations F.12 and F.12), this becomes 

r/, = z {(df.Pf I ^/e"^'^-"^ I C,Pc )( C,Pc\ J AO ■ md'i I ^.,^.,7. ) + 

C,Pc 
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(G.2) 



The first and seventh terms above cancel each other, as do the third and fifth. 
Therefore, we only need to consider terms with energy denominators. We replace 
A by A^^') + A*^^) , where v4*^^) and A'-^'> arc the non-gradient terms in the multipole 
expansion of equation (E.30): 
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{a,,;,D^l(0, - al j-l)^-'D^^l^{0, -^)} , (G.3) 



{ak,xD^iil{0, -d, -v) - al^{-lf-^Dfl^{% -^)} . (G.4) 
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An explicit form for J is also needed, so we write J = J^'^^ + J^p^ , where 

-^^^^ = ^J2{^M-^j),Pj}- (G-6) 



Here fij is the magnetic moment and pj the momentum operator for the j'th 

'fi 



nucleon. To keep track of the sixteen different terms in Ti'^, we write 
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The superscript a2, for example, indicates that this is the term with ■ A'-^\ 
As a reminder, l> = {e/h)^{xp). 

We start by examining the integral / J^°') (^) • A'^^ (O'^'^C- Including only the 
^-dependent terms for now, this integral becomes 



j(-)(^-).i'(l)(^-)rf3^ 



3=n,p 
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X j^{uj^)LYLM{C) ■ lJ-pSp5{S. ~ Xp) + finanS{^ - x„) 



(G.13) 



After performing the integrals, we substitute f = Xp — Xn and neglect recoil 
terms to get 



jt-)(|) .i'(l)(|)rf3^ 
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(G.14) 



Usi ng equ ation ( E.29 ), the spherical Bessel recursion relations (equations K.ll 
and K.12 ) , and the curl version of the gradient formula (equations K.14 - K!T^ 
equation ( G.14 ) becomes 
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{tip - (-l)^A'n) S+{tlp+ (-l)^Mn) t , (G.16) 



where we have defined S = {dp + (T„)/2 and t = {dp — dn)/2. 

We are now ready to calculate the Tj^''^^ term. Removing the center-of-mass 
moti on immediately and inserting equations (F\2), (G^), (G^), and ( G.16| ) into 
(^), we get 
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The t term gives zero contribution since this operator changes the spin, but 
Si = Sc = 'i-- The L + 1 term in the fourth hne wih also be dropped since it is 
of smaller order than the L—1 term, and the L — 1 term turns out to be small 
itself. We now use the fact that 



fjLM{f)-V = [Yl'S)V]j^, 

where V is any vector operator, to get 



(G.18) 
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Again, iPl{x) = + x-^ji^{x). The M sums can be removed to produce the 
final result: 
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The other al terms are very similar so we just write the results here: 
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We turn now to the a2 terms, following the same steps used in the derivation 
of equation (G.16): 
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Putting this into equation (|G.8|) gives (again neglecting the t terms): 
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Using equation (G.18) to convert the dot products into tensor products, we see 
that this becomes 
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The final expression for A4^}-''^^ after removing the M sums is 
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The expressions for the other a2 terms are: 
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Next, we will look at the pi terms. Inserting equations ( |G.3| ) and ( p.6| ) into the 
^-integral gives 
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A second term, which would involve the gradient of the wavefunction, gives no 
contribution and has been dropped. The gradient in the above equation can be 
evaluated with the help of equation ( K.13| ). We find 
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We are now ready to calculate A4^j\'^^. Including only the main {L — 1) term 



from above, we insert two complete sets of f states into equation (p.q) to get 
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After some simplification, this becomes 
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The final expression witliout tlie M sums is 
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We simply list the other pi terms: 
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I L Jc \ f Jc L' 1 

-M/ Ai -Xi + Mf J \ Xi - Mf Mf-Mi-Xi Mi 

{l'll\\ [Yl^i^L]^ \\LcScJc){LcScJc\\Yl' \\in) x 



OO /"OO 



dr dr' r r' u,'{r) Ui{r')tl)^'(^^^)GLc{r,r' ■,EQ)iL_^(^). 



Lastly, we calculate the p2 terms. Again, we simplify the ^-dependent terms 
first: 



jlp)(f) .1(2) (1)^3^ 



(G.41) 



E j {<^Ai~Sj).Pj]-^j^i^oyLM{i)d^^ (G.42) 



S {^'- ■ ["^-(^^^^Mif)] + 2jz,( y)^^*^!^)^. V,} . (G.43) 



Using the gradient formula (equation K.13) and the form of the deuteron wave- 
fimction (equation A.9), we can show that 



and 



Therefore, 



r3.(^)YLM{f) 



V(f I d. 



3..(^f) + ^.;(^f^) Y,.Ar) (G.44) 



E 

I 



r±^y-]ir). 
dr r 



(G.45) 



Jlp)(^-).l(2)(|)rf3^ 



(G.46) 



^Y,M{r) 3,.(-) + -,!(-) + 2,.(-).- ^ , (G.47) 
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where the derivative in the last term acts on a wavefunction. Inserting this into 
equation (G.8), we get 



LcSc Jc A/c 



21 



2md 



X (G.48) 



YL'M'ir) \r){r\C) 



dr' 1-'^ 9 

( C I f) X 



( ^1 -1 



L{L + 1) 2mp 



2i 



This simphfies to 



E EEE 

LcScJcMc W LML'M 



2nip LOf 



flLJi 



2md 



.-«<.,,(o,-,-.,z,-uo,o,o)/i|±||p:= 

( Z'llM/ I Yl'M' I LcScJcMc ){LcScJcMc \ Ylm \ HIM, 



oo />oo 



"'0 



/ dr' r dr u,{r')^L'{'^)G,Jr,r';Eo) x 



After removing the Af sums, this becomes 



CO oo 1+L Jc + Sc 2 

Hr=E E EE E E E 



l=0,2l'=a:2 L=l L' = l Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ 



2mp uif 



(G.50) 



2md 



(2L + 1)(2L' + 1) / 1 



L' 



Jc 



LL'{L + 1){L' + 1) \ ~M f Mf- M, - \ M, + A, 



Jc L 1 



Yl. \\LcScJc){LcScJc\\ Yl ) x 



oc />oo 



"'0 



r' dr' r dr uAr')^L'(^)G,^{r,r';E^) x 



j.GY)n.[r) + ^f^i^)u,{r) + 2rj,{^)u[{r) 
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The other p2 terms are: 



oo oo 



1+L Jc+Sc 2 

7re 



2m„ 



;=0,2 i'=0,2 L=l 1,'=! Sc=0,l Jc=|l-L| Z,c=l Jc-Sc| ^ 



(G.51) 



2md 



:L'+L( i\L' + Jc-\f-\i-Mf-Mi p){L') ^ 



(2L + 1)(2L' + 1) 



1 



Jc 



LU{L + 1){L' + 1) \-Mf Mf-Mi-Xi Mi + Xi 
Jc L 1 



-Mi - Xi Xi Mi 

{LcScJcW Yl liai) 



( mil Yl' \\LcScJc ) X 



oc /■•oo 



/O ^0 



r' dr' r dr Ui{r)ipL{^)Gi^^{r,r'; Eq) x 



OO oc 



1 + L 



^^t^^ = - E E E E E E 'e' ^^^x 



2mr, oji 

l=0,2l'=0,2 L=l L' = l Sc=0,lJc = \l-L\ Lc = \Jc-Sc\ ^ 



Mf-Mi-Xi,-Xf 



(0, 



(2L+1)(2L' + 1) / 1 L 



Jc 



LL'{L + 1){L' + 1) \-Mf Xi Mf-Xi 

Jc L' 1 

-Mf + Xi Mf-M,- Xi Mi 



Yl \\LcScJc ) X 

noo 
r' dr' r dr w,-WV'.(^)G.^(r,r';i;^) x 



oo oo 



1+L 



Jc + Sc 2 

e -IT 



2m„ 



^r=E EEE E E E _ 

/=0,2 /'=0,2 L=l L'=l Sc=0,l Jc = |l-L| Lc = l^c-Sc| 



hiOi 



2md 



i^'+^(_l)^' + ^o-A,-A.^C^O_^^_^^_^^(0,-^,-^)x 



(2X + 1) (2X^ + 1) / 1 L Jc 

LL'{L + 1){L' + 1)\-Mf Xi Mf-Xi 

Jc L' 1 

-Mf + X, Mf- M, - X, M, 



l'n\\ Yl \\LcScJc ) X 
( LcScJcW Yl' IKll r' dr' r dr u,(r')Vi.'(^)G^^(r,r';£;^) x 
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Finally, there are terms in 7^.^ where the substitution A = V<i> for the A in 
the Hamiltonian is not made at all. We can instead substitute one of the other 
two terms in the multipolc expansion. These terms are generally small, so only 
the ones with non-negligible amplitudes are listed below. To keep track of the 
various terms, the names will indicate which combinations of J and A were used 
in its derivation. For example, Ad^J^''^^'" is the term derived from 



E 



( rf, I / .P^i'' ■ A^-' I C )( C I / ,7''"' ■ .I'll I C ) 



-n^f-^-^+E,,-Ec + i 



(G.52) 



IS 



The most important term below is the M 



rl.crl 



term; this is the double Ml 



interaction which dominates the photodisintegration cross-section at threshold. 
Here are the results: 



l+L 



« ^(7l,(Tl,uncr 



'jTe^{?iuif){hu!i) 



2ml 



X (G.53) 



E EEE E E E • 

1=0,2 ;'=0,2Z,=1 L'=l Sc=0,l Jc=\\-L\ Lc=\Jc-Sc\ 



Jc L 1 
-M, - A, A,- Mi 



1 L' Jc 

-Mf Mf-M,-X, M, + A, 

{( ni|| [mp - i-lffin] [Yl'-i ® S]^, WLcScJc ) 
( LcScJcW [fip - (-l)'^Mn] [Yl-i ® S]^ 11/11 ) + 

( z'uil + {-^f^in\ [Yl'-i 1]^, \\LcScJc ) 

( LcScJcW [Mp + (-l)^Mn] [Yl-1 i]i 11/11 )} X 

rdrr' dr' u,{r')3,--.(^)G,^{r,r';E^)3,_,(^)u,{r), 

00 oc l+L ■Tc+Sc \ 

\p \p \p \p 7re (nujf}{noji) 

2-^ .Z^ /l^ 2m^ 

i=0,2 ('=0,2 L = l L' = l Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ P 

»'' + '^'(-l)'^' + ^^-'*"^^A,A/I?i^p_^^_;,^^_;,/O,-0,-0)v/(i+l)(i'+l)X 

/ 1 L Jc W i' 1 

V -M/ Ai Mf-\i Mf Mf-Mi- Xi Mi 

{{ [np - (-l)^Mn] [Yl-1 S]^ WLcScJc ) 

( ic^c Jell [mp - (-l)^Vn] [Yl'-i 5]^, 11/11 ) + 
( /'nil [mp + (-l)^Mn] [Yl-i llic^c Jc ) 

( Lc^cJcll [mp + (-l)"^ Mn] 0^]^, 11/11 )} X 



(G.54) 
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r dr r' dr' u,{T')j,,_,(^)G,^{ry;K)3,.,(^)u,{r), 



10 Jo 



oo oo 



1+L Jc + Sc 



^y^a'2,(j2,uncr 



2ml 



X (G.55) 



E EEE E E E 

1=0,2 I' =0,2 L=l L' = l Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ 
^L+L'^_^f + Jc-M,-M,-X,-Xf^iL')_^_^^_^^^^^ 

1 L' Jc \ f Jc LI 



-Mf Mf-Mi-Xi Mi + Xi J \ -Mi - Xi Xi Mi 
{( mil [np + (-l)^Vn] [Yl' ® S]^, WLcScJc ) 
( LcScJcW [Mp + (-l)^Mn] [Yl ® S]^ \\in ) + 

{ mil - (-l)^'/in] [11. ® t]^, IILc^C ) 

( LcScJcW [^^p - (-i)^Mn] [Yl (s> t]^ win )} X 



DC /■•oo 



JO 



rdrr' dr' u,{r')j,,{^)G,^{r,r';Eo)j,{'^)u,{r), 



OO OO 



M 



cr2,o-2,cr _ 



1+L Jc + Sc 



= EEEEE E E 



\nu}f){nuji) 

2ml 



1=0,2 l'=0,2 L=l !,'=! Sc=0,l Jc=|l-i| ^0=1^0-501 

.L+z,'(_i)Z,'+7o-A.-A,^(i;0_^^_^^ (0, -6, -cl>W{2L + 1){2L' + 1) X 

( I L Jc \( Jc L' 1 

V -Mf A, - A, j 1^ A, - M/ Mf-Mi- Xi Mi 

{{ mil [fip + (-l)^/x„] [Fl ® 5]^ WLcScJc ) 

( Lc^cJcll [mp + (-l)^Vn] [i^L' ® S]^, 11/11 ) + 

( Z'llll [^^p - (-1)^A^„] [Fl (g> t]^ WLcScJc ) 

( LcScJc[[ [mp - (-l)'^Vn] [ii' ® 11/11 )} X 



OO pOO 



r dr r' dr' u,{r')j,,(^)G,^{ry ■,E',)j,(^)u,{r), 



JO 



(G.56) 



M 



pl,2?l,uncr 



1+L Jc+Sc 



= E E E E E E Ex (g-^t) 

i=0,2Z'=0,2 L=l = l Sc=0,l .7c = |l-L| Lc = |Jc-S'c| 
'Ke'^{nUf){nWi)XiXf ^ .L'+Jc-M,-Mi-X.-Xi ^ 

1 L' Jc \( Jc i M 

M/ Mf-M,-X, Mi + Xi J \ -Mi - A, A, M, J 

{ l'n\\ [Yl,., ® L]^,^, WLcScJc ){ LcScJc[[ [Yl-i ® L]^^ ||/11 ) x 
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r dr r' dr' uAr')j,>.,{^)G,^ir,r';Eo)jr,-.i'^)u,ir), 



10 Jo 



oo oo 



1+L Jc + Sc 

Mff'-' = EEEEEE E- 

i=0,2i'=0,2L=lL' = lSc=0,l Jc = \l-L\ Lc = \Jc-Sc\ 



2ml 



^L+L'^_^^L'+Jc-X.-Xf^iL')_^^_^^ _^^^^ -0, -0) [{L + 1){L' + 1)]-* X 

( ^ L Jc \ ( Jc L' I \ 

\ -Mf Xi Mf-Xi J\Xi- Mf Mf-Mi-Xi Mi)'' 

( mil [Yl-1 ® L]^M WLcScJc ){ LcScJcW [Yl'-i ® L]^,m' II^H ) x 
£ £rdr r' dr' u,{r')j,,.,{'^)G,^{r,r';E',)j,_,{^)u,{r), (G.58) 

_ (G.59) 

= \Jc-Sc\ 

j^L+L' f^_'^^L'+Jc-Mf-Mi-Xi-Xf ^ 



oo oo 1+L Jc-^Sc 

Mfr^'— =T.T.T.T. E E >< 

1=0,2 ('=0,2 L = l L' = l Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ 

Tre^(J%ujf){huJi) 



2m^ 



D 



Mf-Mi-Xi,-Xf(^' ^' 



/ (2£ + l)(2L' + l) 
LL'{L + 1){L' + 1) 



1 Jc 
-Mf Mf - M,- X, M, + X 

( /'11|| Fl' WLcScJc ){LcScJc\\ Yl \\in ) X 



Jc L 1 

Mi - Xi Xi Mi 



^0 



dr r' dr' 



2r'jA'^WAr') 



G,^ir,r';Eo)[j.{'^)u,{r)+ 



OO oo 



1+L Jc-\-Sc 



Mf-''- =T.T.T.T.T. E E 

i=0,2;'=0,2 Z,=l L'=l Sc=0,l Jc=|l-i| ic=| Jc-Sc| 



TTe'^{hujf){Uu)i 
2^^^| 



(10.60) 



/ 1 L Jc W i' 1 

-Mf X, Mf - X, ) \ Xi- Mf Mf-M,- X, M, 

{ l'll\\ Yl \\LcScJc ){ LcScJcW 1l' ||/11 ) x 



(2L + 1)(2L' + 1) 



oo poo 



JO 



r dr r dr 



'M^)u.ir') + ^j'A^)u^ir')+ 



'ir'h'(^)<{r') 



G,Ar,r';Eo)[h{^)u,{r)+ 
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(Tl,(T2,uncr 



Jc^Sc 

EEEE E E Ex 

1=0,2 1'=0,2 L=l L'=l Sc=0,l Jc=\l-L\ Lc=\Jc-Sc\ 



2ml 



(-1) 



L'+Jc-Mf-Mi-Xi-\f 



\fD^M^_M,-x,,-x, (0, -0, -ct>W{2L+l){L' + l) X 

1 L' Jc 

Mf Mf - M,- Xi Mi + X 

- (-l)'''Mnl [Yl'-i (S> S]^, WLcScJc ) 



Jc L 1 

Mi - Xi X, Mi 



( LcScJcW [mp + (-l)''Mn] [Yl ^ S]^ lini ) + 
( mil [mp + {-lf^^n] [Yl'-i ® t]^, WLcScJc ) 

{ LcScJcW [/^p - ^ln] [YL(^t]^ Ij/ll )} > 



cx; /■•■DO 



M 



^0 
(T2,crl,uncr 



r dr / dr' u,{r')3,,_,(^)G,^{r,r';Eo)j,{^)u^{r), 



■DC' OO 



1+L Jc+5^c 

E EEE E E Ex 

i=0,2 ;'=0,2 L=l t'=l Sc=0,l Jc=|l-i| Lc = \Jc~Sc\ 

me'^{nujf){nuji) .L'+j^_M,-M,-\,-Xf 



XiD 



(0, -0, -(t>)^y{2L' + l){L + l) X 



{(nil 



2ml 

{L') 

Mf-Mi-\i,-Xf 
1 L' Jc 

Mf Mf - M,- X, M, + X. 

L' 



Jc L 1 

-M, - Xi Xi Mi 



/ip + (-1)^ [Yl' (g> S]^, WLcScJc ) 
( LcScJcW [mp - (-l)'^A*r,] [Yl-1 » S]^ Will ) + 
( riljl [fxp - i-lffin] [Yl' O t]^, WLcScJc > 

( LcScJcW [Mp + i-lff^n] [Yl-1 1]^ \\lll )} x 



OO pOO 



^0 



r dr r' dr' u,{r')j,,{'^)G,Jr,r';Eo)j,.,{^)u,{r), 



OO OO 



M 



(Tl,pl,uncr 



1+L Jc+Sc 
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(G.61) 



(G.62) 



= E EEE E E E 

i=0,2i'=0,2L=l L' = l Sc=0,l Jc = |l-i| ic = l'^c-Sc 



^-{huJf)(hLOi 

2m2 



ji+i''(_l)i'+Jc-M^-M.-A,-A^;^.;^^£,(i') 



1 L' Jc 

-M/ Mf-Mi- Xi Mi + Xi 



M,-M.-A„-A,(0,-^,-'/')y|^X 

Mi - Xi Xi Mi 
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( Z'llll [mp - (-l)^Vn] [Yl'-i ® S]^, WLcScJc ) 
{LcScJcW [Yl-i(^L]^ \\in)x 



-'O 



r dr r' dr' u,{r')j.'-.(^)G,^{r,r' ■,Eo)h-.(^)u^{r), (G.63) 



A^pl,a'l,uncr 



aL+L' 



CO OO 1+1/ Jc'-\-Sc 

E EEE E E E 

i=0,2i'=0,2L=l L' = l Sc=0,l Jc = |l-i| ic = |Jc-5'c| 



2ml 



L + 1 
L' + l 



Jc L 1 
Mi - Xi Xi Mi 



1 L' Jc 

-Mf Mf - M,- Xi Mi + X 

{l'n\\ [YL'-i(g>L]^, WLcScJc) 

( LcScJcW [mp - (-l)''Mn] [Yl-1 « S]^ Win ) X 



OO fCC 



JO 



rdrr' dr' uAr')].>-^(^)G,^^{r,r' ■EQ)3,_,(^)u,{r), (G.64) 



OC' OO 



2ml 



2L + 1 
L'+l 



l+L Jc+Sc 

^pi,.2,.nc. = j2j:j:j:j: e e 

1=0,2 i'=0,2 L=l L'=l Sc=0,l Jc=\l-L\ Lc=\Jc-Sc\ 

[ 1 L' Jc \( Jc LI 

\ -Mf Mf-Mi-Xi Mi + Xi J \ -Mi - Xi Xi Mi 

{l'n\\ [Yl'-i®L]^, WLcScJc) 

{ LcScJcW [Mp + (-l)''Mn] [Yl ® S]^ IIZll ) x 

rdr r' dr' uAr')j.'-.{'^)G,^{ry;Eo)j4'^)u,{r), (G.65) 



(72,pl,uncr 



l+L Jc+Sc 



EEEE E 

l=0,2l'=0,2L=l L'=l Sc=0,l Jc=|l-i| ic=|Jc-Sc| 



2L' + 1 



-Mf Mf-M,-X, M, + X, 



Jc L 1 

-M,-X, A, M, 



{ mil [/ip + (-1)'' Mn] ® S]^, WLcScJc > 
(Lc^cJcll [Yl-i(E>L]^ lini) X 



/•OO /-DC / 

^c!rr'dr't.a^')i.'(^)G.^(r,r';£;o)j.-.(^K(r). (G.66) 



Appendix H 

Recoil Corrections 



We will now return to the recoil corrections to the transition matrices which arise 
from the P operators in equations ( F.12| ) and (F.13). The main contributions 
from these terms are 



T 



cm.uncr 



(H.l) 



E 

C,Pc 



E 

C,Pc 



df,Pf\ ^V<I>/(f)e-'=/-« - F I C\Pc){ C,Pc 



Jki-R 



di , Pi 



hw, + Ed, -Ec- Pll'2.mi + ie 



df.Pf I ^V<i>,(f)e*-^-P| C,Pc)(C,Fc 



(H.2) 

e-^^V-^ I d,,P, ) 



-TiLjf + Ed, - Ec - Pcf^md + ie 



Since we are working in the lab frame, P \ Pi ) = 0. After integrating out the 

-CIT 



center-of-mass variables and using A — V$, 7''™'""'=' becomes 



cm.uncr 



2Trfi 5{P, + h- Pf - k 



J J 



V Vujf V 



(H.3) 



E 



C - ^ ^d, - Ec + ie 

The commutator can now be evaluated; this produces two terms: 



t: 



cm.uncr 



2Trh 6{P, + ki~ Pf- kf) 



where 



E 

c 



V v^f 



,T-cm,uncrl . rT-cm,uncr2 



,(H.4) 



flLOi 



2md 



TlOJi 



2mj 



Ed, - Ec + ie 
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flip 



ieh , 



(H.5) 



Starting with 7^';™>""'='i^ -^^g insert equation (F.2) for <l>i and expand the expo- 
nential into partial waves to get 



oo oo 



MT--' = E E E E X 

C L = l L'=0 MM' 

df I '-^{h ■ el^)M-^fj.'i'^)Yl,M'Ckf)YL'M'if) I C : 



(H.6) 



1 



{C\ 



TlLJi 



2m.d 



4fl (0,0,0) 



z^+i 27r(2L + l) u^r 

T^V L{L + 1) M — )YLM{r)\dn, 



where the T^^ has been replaced w ith -M fj, according to equation (A_/7). The 
function t/i^ is defined in equation (F.25). Next, we insert two complete sets of 
radial states: 



oo oo 



* ^cm,uncrl 



-E E EEE(^-^A,) 

W LcScJcMc L=l L'=0 MM' 



hujj 



into,)' 



2md 



(H<7) 



{ I'UMf I Yl'm' I LcScJcMc )( LcScJcMc \ Ylm \ HIM, ) x 



/27r(2L+ 1) 



OO /'OO 



r dr r' dr' uAr')j,'(^)G,^{ry;E^)^,(^)u,{r). 
The final answer is then obtained by removing the M sums : 



M 



l+L JcVSc oo oo 

E E E E E EEx 

(=0,2 l'=0,2 Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ L=l L'=0 



47re^ 



[flUJi 



2md 



(-1) 



Jc-Mf-M.-X, 



,l^l'y* , 27r(2£+l) 

( Fl' WLcScJc ){LcScJc\\ Yl ) x 



1 L' Jc 

-Mf Mf - M, - A, M, + A, 



Jc 



L 1 



/•OO /'OC 



(H.8) 



r dr r' dr' u,ir')j,,i^)G,Jr,r';Eo)M'^)Mr)- 
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We just list the results for the other MJf terms: 

L+L' oo oo 



ji^cm,uncr2 



M 



cm.crl 



_^cm,cr2 



(H.9) 



i=0,2 l'=0,2 L=\L-L'\ L=l L'=0 



•L — L'-w-* 



(2i' + l)(2L + l) 



1) 



( mil 11/11 



L V L 



L 



L' 



Q J \ \i Mf-Mi-Xi Mi-Mf 



-Mf Mf - Mi M., 

l+L Jc+Sc oo oo 

E E E E E EEx 

(=0,2 i'=0,2Sc=0,l Jc = \l-L\ Lc = \Jc-Sc\ L' = l L=0 



(H.IO) 



TujOf + 



2md 



(^_l^L'-\f-Jc 



/ (2L + 1)(2L- + 1) (,.0 . 
( mil Yl WLcScJc )( LcScJcW Yl' ||/11 ) x 



1 L Jc 

-Mf Mf 



Jc L' 1 

-M/ Mf - Mi Mi 



nOO / 
r dr r' dr' u,{r')i,,,(^)G,^{r,r' ■,E',)j,{^)u,{r), 



L+L' 



E E EEE(^/-aJ^(-i)^'— ^ 

L=\L-L'\ '=0,2 Z'=0,2 Z,=0 L'=l 



Tod 



(H.11) 



i^+^'(2L' + 1)(2L + 1)^ ^|±i^<)_,^ (0, -<i>) X 
<^'ll|l^^ll^l^)(o o)(o Mf'-M. M^Mf 



1 L 

-M/ Mf - Mi M, 



Appendix I 

Relativistic Corrections 



The relativistic correction to the charge density is given by 



P 



eh 



4m^ 



(I.l) 



where k is the anomalous magnetic moment. We can calculate the effect of this 
correction by inserting 

into equation ( |F.9| ). Using the approximation that V$i w A = e\^e^^^'^ , the 
integrals over the dummy variables can be performed immediately: 



(1.3) 
(1.4) 

(1.5) 
(1.6) 



/J' 



where 



{2Kp + l)e*^--^-/2 - 2K„e-*^'-*'/2 
{2Kp + l)e-^'''ff/^ _ 2Kr^e''''■^/^ 



CA, ■ (p X S), 



and S* = ((Tp + 0Vi)/2. The term proportional to i = (iTp — o>i)/2 is zero if we 
only work to first order in the relativistic correction, and therefore has not been 
included in the definition of ^. To this order, recoil corrections can also be 
neglected. Thus, the correction to the scattering amplitude is given by 



(1.7) 
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where 



^rell _ ^ {df\[H-P,^f]\C){C\[H"P,h]\di) 



c 



M}f = --{df\^f[H^P,h]\di) + -{df\[H-P,^f]h\d.)- (1.13) 



j^rel2 _ ^ {df\[H-P,^^f]\C){C\[H-P,^i]\di) 

- ^ nco, + E,,-Ec + ie ' ^'■^> 

^relS _ ^ {df\[H-P,^,]\C){C\[H-P,^f]\d,) 

= -ru^f + Ea,-Ec+ie ' (^-'^^ 

We now evaluate the commutators in the above terms, using, for example, the 
fact that 

( df I [H-P, ^f]\C) = {df\ (Ea^ - Ec)^f I C ) (1.12) 
to cancel some of the energy denominators. This results in 

X/ , , i * , , N -.^sr^ {df\^f\C){C \ iii\d,) 

c 

M}f = -l{df\ ^f[H-P, h] \di) + l{df\ [H-P, \di)- (1.14) 

X, , , ^ i , , ^ n2V- ( '^f I *f I C'li C I I ) 

in..)i d, I I d. ) + (/.u..)^ E Z/+^,:-Ec + ie ' 
A^/f = —{dfl */] I rf, ) + ^( rf/ I [^f"^ ^i]^f \di)+ (1.15) 

A^/f = -l{df\h[H^',^f]\di) + l{df\[H^P,h]^f\di)+ (1.16) 
(;ia;/)( I I ) + (;ia;/) 

The largest terms are the third terms of each group. We call these the "contact" 
terms and sum them to get 

^rei.cont ^ _ ^^.^^ | ^.^^ ^ | ^ . ^ ^ ^^^^^^ | j^.^^^j | ^ . ^ _ 

(na;i)(d/| (L17) 

The first term is Lu/rrid times smaller than the others, so we will neglect it for 
now. The commutator in the dominant terms is basically 

[p, = -eiV^i « -eiex.e''^'-''/^, (1.18) 
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so that these terms are equal to 

2Kn{hLjfe~'^-^/^ -huj.e'^-^/^)] \d,), (1.19) 



where q=kf — ki and K = {kf + ki)/2. Note that the correction which depends 
on Kp is of order uj /nid times larger than the k„ piece. 
This term is easily evaluated; the result is 



2 2 oo 

= -^E E E EEx 

P i=0,2 i'=0,2i=-l,04 J=o -f-=o 

/ L 1 J A / 1 J I \ 

\ -j + Mf ~ M, j M, - Mf J V Mf - M, M, J 

t^(-iy-^'^YlM{qW2J+l( rilll [Yl ® S]j Will ). (1.20) 



We now examine the terms in equations (I.13)-(|3^) which contain commuta- 
tors. Adding these eight terms together gives 

■' ZrUp 

<->{<l^,#,}. (1.21) 

Since we are working in the transverse gauge where V • A = 0, the V^'I' terms 
can be neglected . Also, looking back at the definitions of we notice that 

^2 

V'^-.m = -J^^^(f), (1-22) 
Therefore, the commutator terms become 

(1.23) 

This is a correction to the contact term. We can account for the "non-commutator" 
part of the above equation by making the substitutions 

hco, ^ hco,-^^^ (1.24) 
Zirip 

hcof - ^f+^^^ (1-25) 
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in equation ( |I.20| ). This amounts to a correction of order (aj/m^)^ smaller than 
the main term, so it will be discarded. We add the (previously neglected) first 
term of equation ( |1.17 ) to the remaining part of the commutator to get 



M 



re/,comm 



hujf — huji 



df I ^j^-/ + ^/^-i I 



(1.26) 



This is the next-to-leading-order term, along with the k„ piece. It is assumed 
to be small at energies below 100 MeV and is not evaluated further. 



Appendix J 

Pion Terms 



We now consider diagrams, such as those in Figure J.l, where a pion is ex- 
changed between the two ■^N vertices. The first terms to be calculated are 
those which are needed to satisfy the low-energy theorem described in Section 
3.6. We then calculate the potential energy portion of the double commutator 
term in equation (F.21); this exactly cancels the other pion-exchange diagrams 
at low energy. Lastly, a correction to the vertex is presented. 

The Hamiltonian for a point interaction between a pion, photon, and nuclcon 

is 



E 



(J.l) 



where ^ « 0.075 and fj is the isospin operator for the jth nucleon (the tilde 
indicates a vector in isospin space). A{xj) is given by equation ( |A.2D and a 
similar expansion for the pion field is 



q 



-iq-Xj 



(J.2) 



where a]^ - creates a tt^ with momentum q. Inserting the above equations into 
equation (F\l) gives 

(J.3) 



rriT /-7-7rl,uncr , /-T^7r2,uncr . /t— Trl.cr , /T^7r2,cr 



where 



EE 



C,Pc 



V(fq pe^ 2tt 2tt 

(2^)3 ml V^^uJJUr.E^ 



(J.4) 



df,Pf\{<f2-ei;, 



-, — ik 



fX2_ 



V2 



\C,Pc ) X 



hoji - -B^ + Ed, - Ec~ Pl:/2mD + ie 
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{C,Pc\ • EAje^^" -~ii=ri^e-'«-i \ d,,P, ), 

^.2,uncr _ ^ ^ / Vd^q pe^ 2^ 2^ 

( df,Pf I (ai . e*^)e-^~/-^^i=ri±e^^-~^ | C, Pc ) X 
1 



- -Bvr + - Ec - P^l2mD + 
1 



- £^,r + Ed, - Ec- Pl,/2mD + ie 
( C,Pc| (^1 • 6%)e-^'/-^^-^ri^e-«-~^ jd^P, ), 
^.2,cr ^ /■ Vd^q Pe^ 2^ 2^ 

C,Pc 

1 



-Twof - + Edi - Ec - P^/2mD + is 

( C,Pc| (^2 • e%)e-^^/-^^-^r2^e-'--~^ |d„P ). 

These four terms correspond to the four time-ordered Feynman diagrams of 
Figure [J.l[ . We begin with 7^'^^^""'^^ After removing the center-of-mass states 

from 'Y'"^''^'^" and inserting complete sets of r and f states, we get 



, d\ npe^ 
(27r)3 ml E^ 



_^^l,u„cr 



yy I d\ d\'^ 



r ) X 



{df\ (a2-el,)e^'^^-'-r2±e 

^'^\Lo,^E,+Ed,-Ec -Pi /2mD + ze ^ ^ ' ^ 
( f| (ai • eAje*"'-Vi^e-^«-^^ |d, ). (J.8) 



We now make the assumption that 



E^ = J ml + » huj, + Ed, - Ec- Pll2mu. (J.9) 
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This collapses the sum over C and the integrals over r' so that 



(2^)3 ml + 

Performing the three-dimensional integral over q yields: 

( ^1 (^2 • ^*Xf)i^l ■ {T2-TI+ + T2+Ti_) C^^'^" ), 

where K = kf + ki. 

The isospin operators are evaluated next. We can write 

r2-ri+ + r2+ri_ = 2 (f i • ^2 + 1) (J.12) 
= 4r2-Tf-T| + 2 (J. 13) 

= 4r2-4, (J. 14) 

where we have defined the total isospin operator T = (fi + T2)/2. Since the 
deuteron has zero isospin, Mf^ becomes 

<7rl,uncr _ / ^j3^J / J I ^\ / ^1 f* ; 



J TTijj. r 

(J- 15) 

Inserting the partial wave expansion for the exponential (equation A.lCj ) and 
rewriting the dot products as sums over spherical components yields 

A^;-""" = E E E E E Ex 

L=0 M=-L 1=0,2 ;'=0,2 j = -l,0,l j' = -l,0,l 

-J^^\-ly+= Y£^{K) / dru4r)^—j,{ — )u,{r) x 

(eA,)-.(eAj/( /'llAf/ I a2,,(7i,_,vyiM I lllM, ). 

We evaluate the general matrix element 

( L.S.J^Ah I <y2,j<Ji,-rYi„,Yi,^, I L2S2J2M2 ), (J.17) 

since this form appears in all of the pion terms to follow. The first step is the 
formation of tensor products using equation (K.l): 

a2,,ai,_,. = ^ (-1)^^' V25^ ( J ) (^2 ® ai)^,^, (J.18) 



S'M' 
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and 

JM S'M' LM 

I I' L \ f I I' L \ f I 1 S' \ f S' L J 
j i TO to' -M ) [j -f -M' ) [ M' M ~M 



/ (2; + l)(2r + 1)(2L + 1)(2S" + 1)(2J + 1) . , 
Y ^ 4^ [(^2 ® cr,)s, ® Yl\jm ■ 

The sums run over all possible values allowed by the properties of the Sj symbols. 
Next, we use the Wigner-Eckart theorem. This gives 

( L^SiJiMi I a2^jai^-rYirnYi,rn' I L2S2J2M2 ) (J.20) 



JM S'M' LM 

I V L \ f S' L J \ f Ji J J: 



1 1 S' \( I V L 
J -f -M' j I 



2 

TO to' -M j \ M' M -M J V -Ml M M2 

1 (21 + 1)(2;' + 1)(2£ + l)(2y + 1)(2J + 1) 
V 47r ^ 

( Li^i Jill [(ct2 ® cri)s, ® FlIj 11^25*2 J2 ). 

This tensor product can be separated into spin and orbital angular momentum 
parts using equation (K.3). Thus, the matrix element becomes 

( LiS'i JiMi I Cr2jO'l -j-'YimFi'm' I L2S2J2M2 ) 

= E E ^.{-1)^''+''-'^+^'+'^-''^+'' X (J.21) 

JM S'M' LM 



(2J I 1) / (2^ + l)(2^^ + l)(2^ + l)(2^i + l)(2J2 + l)(2y + l) 
V 47r 
1 1 S' \ f I I' L \ f I I' L \ f S' L J 

j -f -M' ) yo ) ym to' -M ) y M' M -M 

{X i J t |}(5i||(..«.i)„||..)(Li||r,||L.) 

- EEE(-i)''"'^''^'^""'^^''(2^~+i)(25' + i)(o S 0)^ 

/ /' L \f S' L J \ / Ji J J2 
m m! ~M ) \ M' M ~M ) \ -Mi M M2 




{21 + l){2l' + 1)(2L + l)(2 Ji + l)(2J2 + l) 

X 



47r 
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t\\s'){s'\\ S + t\\S2), 



where S = (cti + (T2)/2 and t = (ai — (?2)/2. The orbital angular momentum 
matrix elements can be evaluated using equation (K.5), while the spin matrix 
elements are listed below: 



Si\\s\\s') = VeSs'iSs.i, 

{Si\\t\\s') = VsiSs'oSs^i ~ Ss'i6sio)- 



(J. 23) 
(J.24) 



Using this formula for the specific case at hand, we see that M'^^''^^'^^ becomes 



M 



EEEEEE 

LM W jj' jmS'M' s' 



2^2 



z^(2J+l)(25' + l)(el^)_,(6,J,, / dru,{r) 



— X 




1 J 1 

-Mj M M, 



1'\\Yl\\1){Q6s'i+?>5s'o)- 



The M sums are removed in the final result: 



(J.25) 



M 



E E E E E E Ex 

i=0,2 i'=0,2 j = -l,0,l i=-l,0,l S'=0,2 J=o.l,2 ^=0 



mt 



1 1 S' 

i -/ i' - i 
1 J 1 

-Mf Mf - M, M, 



r 2 



S" i J 

j-y Mf-M,-j+j' M,-Mf 

1 /' 1 
1 I 1 
S' L J 



(J. 26) 



(r||yL||0(-'5s'o + 2(55'2). 



We now turn to Al^^'""'^'^. Comparing equations ( |J.4| ) and (L5), we see the 
only difference is that the nucleon labels 1 and 2 have been switched. Since the 



numbering of the nucleons can't matter, we must have Ai 



M 



fi 



In fact, under the approximation that the pion energy dominates the energy 
denominators, all four time-orderings given by equations (J.4)-(J.7) are equal. 
Therefore, we get 



A^TrljUncr , »^7r2,uncr 



(J.27) 
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Figure 3.2: Feynman Diagram with one ^jmr vertex 



oc »2 2 

E E E E E E E ^^^-d---- x 

i=0,2 i'=0,2 j=-l,04 j = -l,0,l S'=0,2 ,7=0,1.2 -^^=0 ^ 

{2J+1){2S' + l)YlM,-M.-,+fiK) dr u,{r)-- — )u,(r) x 

1 1 S' \ f S' L J 

J -f f -J J\J- f Mf - M, - J + f A-U - M/ 

^ ' 1 I 1 ) {l'\\YL\\l){-5s'o + ^5s'2)■ 
S' L J 



1 J 



Next, we will consider the diagram in Figure J. 2, The Hamiltonian for the 
photon/pion vertex is 



H 



777777 



d^xe'^A^{x) [(t)+{x)(t)-{x) + <j)-{x)<j)+{x)] 



(J. 28) 



and the Hamiltonian for the ttN vertex is 

hf 



(J. 29) 



i=i,2 



The different time-orderings for this diagram are shown in Figure J. 3. Careful 
analysis shows that the sum of Figures J. 3(b) and (c) equals Figure [J.3|(a). There 



are three other time-orderings not pictured which are identical to Figure J. 3 



The net result of this counting is that Figure |Jj(a) needs to be multiplied by a 
factor of 4. Calling the intermediate TV TV state between the two nN vertices C, 



the diagram in Figure J. 3 (a) can be written in third-order perturbation theory 
as: 



C,Pc 



V(Pqi V(Pq2 



1 



(27r)3 (27r)3 Ed^ ~ Ec + huj, -huf ~ P^/2md - E„^ + ie 



df,Pf,lf\ ^ E i'^r^j) (rj-H^j)) |C,Pc,7r2,7/ ) x 
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Figure J. 3: Time-ordered diagrams for Figure 
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C, Pc,7r2,7/ I 



1 



d^x e^A\x) [(j)+{x)(f)^{x)+ 



1 



Ec - P^/2md - E^, + ie 



xr))\d,,P,n^)■ (J.30) 



j'=n,p 



Now, we insert equation (J. 2) for and equation (A. 2) for A. Under the usual 
assumption that the energy denominators are dominated by the pion energies 
(which ehminates the sum over C), 7^^ becomes 

Vd^qi Vd^q2 1 



4^y d^x^ 



(2^)3 (27r)3 -E^ 



-{df,Pf\ 



{(ai • 92 ) e'?--~iri± + (^2 • 92) e^«^-"^V2±} x 



A-kTi V ^J^^^J^i 
1 



o»(9i-92)-2' 



V .JE^. 



-Em VVy rn. 
('?2-'fi)e-'«^-^V2^} \d,,P, 



{{^1 ■ qi) 



-iqi-xi 



(J.31) 



Next we calculate the effect of the isospin operators. Combinations involving 
only one of the nucleons (such as ri+ri_) are already included in the polar- 
izability terms and therefore are neglected here. The remaining terms can be 
evaluated using equation ( J. 14 ); the result after removing the center-of-mass 
states is 



Mf 



d^rd^x£^£Z'-^^^iil,-ix,y 



(27r)3 (27r)3 ml 

( df I r){ r I {(ai • q2){ff2 ■ q'i)e'9--(-72-5)em-(r72+.-)g-*?-5 + ^ 2)} |d, ), 

where q=kf — ki. To simplify the calculation, we shift the variable x by — r/2. 
The integrations over qi and q2 can then easily be performed. Altogether, we 
get 



Mf 



d-^r d^x 



2 — "^Af ■ caJ df r }e 



■iq-{x—r/2) 



(J. 33) 



f\ (ai-Vr) (T2 • (Vr+V^) 



^—ra^xjTi ^ — 7nT^\r—x\/7i 

Anx Anlr — x\ 



The gradient with respect to x can most easily be evaluated by integrating by 
parts. Before evaluating the other gradient, we expand the exponentials into 
partial waves, and also use 



-mTT \r—x\/h 
\f— x\ 



47rm^ 



E E 

L=0 M=-L 



'^)k.(^)YLM{f)Yl^{x\ (J.34) 
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where t'<(>) is the lcsscr(greater) of r and .x, and and arc the spherical 
Bessel functions with imaginary arguments. Also, the integration over the an- 
gular part of X can be carried out because of the orthogonality of the spherical 
harmonics. Thus, A4J^ becomes 

^2^2 ^2 



■' "'^■^ LM L'M' 

{ r\ YlM{q)yi'M'iQ)yLM{r)jr.{Y) J dx x^j^,{qx)iai ■ V^) x 

a2-{Vr +iq) i.'{^^)k,,{'^)- FL'M'(r) \di ). 

L J 71 n X 

The (7 dot products are expanded into spherical components: 

"'"^'^ LM L'M' jj' 

{ df I r )( f| ai,,va2,-,y£M(9)>7'M'('?)>^LM(r)j.(|^) X 

(V_yV, + i9jV-/)ri,M'(r)/L'(9,r) ), 



where 

fL'{q,r)= / rfa; a;V(9a;)«i.'(^^)fci-'(^^) • (J-37) 

Jo ft ft X 

Two applications of the gradient formula are needed to evaluate this expression. 
We write the final result as 



^ A 2 2 £2 oo L oo L' 2 

^f' = -^^(^X-'^^)E E E E E E ExJ-38) 

^ L=0 M=-L L'=0 M' = -L' j=-lfi,l j' = -l,0dn=-2 



Jo 



dr u,{r) u,{r)n{\)ct^M'n'Ot){T)fL'{q,r). 



where 



Ci-'2',y ^ (-l)i-V^'(i' + l)(^^',_y j,)x (J.39) 



L' - 2 L' - 1 1 
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Figure J. 4: Feynman diagram corresponding to T, 



0[V{r) 



L' L' + 1 1 

M'-f+j -M' + j' -J 

d L' + 2\ / d V 
dr r J \dr r 

L' -1 L' l\( L' L' -1 1 

M'-f -M' f)[M'-f+j -M'+f -J 

d L' d L' + l 

dr r J \ dr r 

L' + l L' 1 



(-l)-^L' 



(-l)i-V(i' + l)(i' + 2) 



L' + 2 L' + l 

M'-j'+j -M'+f 



M' - f -M' f 

L' + l L' 1 
M' - f M' f 

1 

-j 



d_ 

dr 
d 
dr 
d 
dr 
d 
dr 





(- 




\dr 


L' + 




r 








v) 


1 


L' + 


1\ 



L' + l 



d 
dr 



r 



(J.42) 
(J.43) 

(J.44) 
(J.45) 
(J.46) 
(J.47) 



We also consider diagrams such as the one shown in Figure J. 4. Just as in the 
previous term, there are 6 time orderings; this contributes a factor of 4 to the 
calculation. The Hamiltonian for the single photon/pion interaction is 

H^^'' = ^ / d'^iqi + 92) • A{x) [0_(gi)0+(q2) - 0+(gi)0_(g2)] , (J.48) 
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which means that this term is 



^fi 2^ J (27r)3 (27r)3 

C,Pc 

1 



- Ec + huJi-?iiOf - P^/2md - E^^ + ie 



( df,Pf,lf\ -f ^ • V,) (fj-m)) \C\Pc,7T2,7f ) 

J=l>2 

( C, Pc, 7T2, lf\£^J d^xiqx + 92) • A{x) [(/._(gi)</.+ (g2)+ 



\C,Pc,.ulf )nu,^ + Ed,-Ec-Py2md-E^^+ie 
{ C,Pc,ni,ji I — ^ (ct„ • i*(f„)) (f„ • \di,Pi,^i ). 



n— 1.2 



This calculation follows the same steps as the previous one. We first simplify 
the energy denominators, and insert the expressions for A and (j): 



50) 



g-i9i.xigifei-xi^^^^ 



{(a^i • 92) e'«-^^Vi± + (a2 • 92) e'«-^^V2±} x 

(^2 • e-'9--"^r2^e*-^} |di,P^ ). 

The only isospin factors which appear are ^_|_ ri±T2:)i = —4. Removing the 
center-of-mass contributions and adding in the crossed term then gives 

= -/d3.d3,AL|g| 16.eW 1 ^\^ id,\f) (J.51) 

{ f\ {(ai • q2){a2 ■ ex.M + Q2) ■ e*^e^«'-(^+^^/2)e^«-(^^/2-^)e-^'=/-^^e-'^^-^^/2+ 

(^2 • 92)(ai • hMQl + Ql) ■ e^^e'9i-(^-^72)ei92-(-r72-x)g-ifc}.Xgifc,.r72 ^ 

(cti • 92)(ct2 • el^Wi + <?2) • f:;,^e'ft-(^+^'/2)e'9-=-(^72-S)eifc;.Xg%.r72 ^ 
(^2 • g2)(ai • el^Wi + <?2) • e;,^e*?i-(---72)ei92-(-r72-x)gifc;.Xg-ife,.72| 

We again shift x by — r/2 for convenience, and rewrite the momentum operators 
as gradients. The 9*1 and 92 integrals can then be performed to yield 

= / d'x—-^{ df I r )x (J.52) 
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Attx Airlr — x\ 



We integrate by parts to evaluate the x gradient, and find it to be zero. Next, we 
expand the exponentials into partial waves and break up the dot products into 
spherical components. After integrating over x to remove two of the spherical 
harmonics, the matrix element becomes 



- / d^: 



oo L oo L 

E E E E E E E 

L=0 M=-L L'=0 M' = -L' k=-1,0,1 m=-1,0,1 i^=-1,0,1 



( r\ rLAf(r)[(-l)^'2y£,jv/,(fc/)fTi,^a2,,.(V^)_^(V,.)KXL'Af'(r)5L'(w/,r) + 



qr ^ 



where 



9L'(uJ,r) = dx X — - — jfcL'i^ — ) — 



-j^,{ujx). 



(J.54) 



h h ' X 

The gradients can be evaluated using the gradient formula to give the final result 

oo L oc L 

Mft--T. i^-^^) 

1=0,2 l'=(),2 L=Q M=-L L'=0 M' = -L' k=-1,04 /j=-1,0,1 i/=-1,0,1 n=-2,0,2 

i^|^^^+^'(-l)''+^+''(^AJ-.(6X.)-.r£M(9)^°°dr u^ir)M^)u,ir) x 

[( I'llMf I {-lf(Jl,^<J2,.YLM{f)YL'+nM'+u-^.{f) \ lUAh ) X 



( I'llMj I ai,^,(72,r.YLM{f)YL'+nM'+u-,.{f) \ HIM, ) X 

Yl,M'{k)Ct^t'.A^h^)9L'i.^^.r)]. 

The C(") and O^") factors are defined by equations (|j.39|) -( |j.47| ), and the an- 
gular element can be evaluated using equation ( J. 22 ). 

We now calculate the diagram pictured in Figure J. 5. There are 12 separate 
time-orderings contributing a total factor of 4, as well as a factor of 2 for the 
term where nucleons 1 and 2 are switched, so the matrix element is 



^7r7 



Vd^qi Vd^q2 Vd^qs 



C.Pc 



|3 (27r)3 (27r)3 



(J. 56) 



( df,Pf,lf I [dj ■ \'j){fj ■ (l){xj)) I C, Pc, 7/, 7^3 ) X 

771-77 
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Figure J. 5: Feynman diagram corresponding to jVl jj. 



Edi + hu>i - Ec - Pc/'^rud - E^^-huf + i 



IS 



±e 



C, Pelf, -^3 I -^iq2 + qs) ■ A(y)(j)^(q2)(j)±(q3) \ C,Pc,tt2) x 



Ed, + hiOi -Ec- P^/2md - E^^ + is 

^ ±e ^ 

( C,Pc,n2 I ^(9*1+92) • A{x)4>^{qi)4>±{q2) \ C,Pc,-Ki,ii ) x 



Ed, -Ec- Pll^rud - E^, + is 

^ f 

{ C,Pc, 7ri,7i I {aj ■ Vj){fj ■ 4>{xj)) \ duPun ). 

UItt 

The energy denominators can be simplified through the assumption that the 
pion energy is the dominant factor. After inserting the expressions for <j) and A 
and removing the overall center-of-mass factors, we get 



r X 



( df I f)( f I { T2,^n,±e"^-^e-"^ry x 



(J.57) 



• gi)(CT2 • qaiqi + 92) • ex,{q2 + qa) ■ e^^ x 

In the calculations of the other pion terms, the next step was to perform the in- 
tegrals over the pion momenta. However, in this case, this turns out to be more 
difficult because of (what will become) the four gradients. We will therefore 
integrate first over the dummy variables x and y. This generates two momen- 
tum delta functions, allowing two of the remaining integrals to be performed 
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immediately. The result of these integrations is 



E 



r{df\r){r\ — T2,^ri,±e /x(J.58) 



[^2 • {q2 - fc/)] [<?i • {q2 - h)] 2(92 • eAj2(gi • el,)e-*«^-*' x 
1 1 1 I , 

The isospin factor is the usual ti^^t\,± = —4. Rewriting the q2 vectors as 
gradients, the matrix element becomes 



(J. 59) 



+ {(fc.,eA.)^(-fc/,el,)}, 



where 



d3<72 



(J. 60) 

and where the gradient only acts on T{ki,kf]r). Defining Z{q,Lo) = (mj + 
h^q^ + Lj"^) / 2h^ quj (as in |52 ) and expanding the exponential, this integral 
becomes 



i-it r ^ ^ f '^'?2 dq'2 ^ -A {~i)^n 

J-.(g2r)riM(r)F£M('?2) 



(J.61) 



[Z(92, t^j) - 92 • fci] [Z{q2,ujf) - q2 ■ kf] [ml + h g|] 



We now use the identity 



00 Q 



= E E 47rQ„W>^a/3(a)r„*^(6), 



(J.62) 



which is true for t > 1, where (3(t) is the irregular solution of Legendre's equa- 
tion. This gives 



m,kf;r) - Y^Y.Y.Y. h'i'^^LMi42)Y:p{<i2)Y:,p,{q2) 



a/3 a'P' LM 

^'?2 ,9 ^2 ■?^ 

Qa[Ziq2,UJi)]Qa' [Z{q2, w/)] 



jL{q2r)YLM{f)Ya0{ki)Ya>f}'{kf) x 

(J.63) 
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The angular integral can be performed by combining two of the spherical har- 
monics and then using the orthogonality relation. The result is 



a(3 a'f3' LM 



a a' L 






a' 


M 




/3' 


M y 



(J.64) 



(2a + l)(2a' + 1)(2L + 1) 2{-iY 



47r 



UiUJf 



■YLM{f)Y^(mY^'P'{kj)Kf^^. 



where 



1 

¥ 



g2 ^^^'^Zl ^ Qa[Z{q2,UJi)]Qa'[Z{q2,UJ f)]. 



(J. 65) 



This notation allows a more compact formulation of the gradient formula: 



V^yiM(r)Ai"i,(a;„c./;r) = ^ (-1) 



p+L+Af 



a=±l 



L + a L 1 



X (J. 66) 



\/i + ^(1 + «)i^L+l,A./+a(r)Ai"+,'L(^., r). 



The final expression is cumbersome and therefore will not be explicitly written. 
Instead, we use eq uatio n ( J.5£ ) and either 2, 3, or 4 applications of the gradient 
formula (equation [J.66 ). Switching the initial and final photons simply causes 
the function T to be multiplied by a factor of (—1)"+" = (—1)-^. 

We would now like to calculate the contribution to the double commutator 
term of equation (F.21) from the pion-exchange potential. This potential can 
be written as 



K(r) = 
where 



1 



1 /e-^ 



47r 



5*12 
X 



3((Tl • f)((T2 • f) - ((Tl • (T2), 



The only operators in the double commutator 



Sir) 
\j.67) 

(J. 68) 
(J. 69) 

(J.70) 



which are non-commuting are the isospin operators. We must therefore explic- 
itly indicate the isospin dependence of l>, which is really a function of Xp. Thus, 
the commutator is rewritten as 

^ E [[(n-f2),r.,,],r.,™]T/— (f)$,(fz)a>y(f„). (J.71) 

/,m— n,p 
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Using the well-known commutator relation 

[Ti,Tj] = 2ieijkTk, (J-72) 

this becomes 

- ^Vr^ir)M^p)^f{^p)- (J-73) 
n 

The term we are interested in is just this commutator sandwiched between the 
initi al an d final deuteron states. By inserting the definitions of $ (equations |F.2| 
and |F.3| ), we find the scattering amplitude: 

L=0 M=-L L' = -AI' M' ;=0.2;'=0,2 ^ ^' 



D,, (0, 0, Q)Dm' - a, (0, -O.-c^y J LL'{L + l)iL' + l) 



where 



Ill,,Acof,u;.) ^ Jdru,{r)uAr){l + r^)j,i^)x (J.75) 

(l + r^)j.K^)5(3 + 3x + x^), 
Ill,u,{u;f,oo.) ^ Jdru,{r)u,{r){l + r-^)U'^)x (J.76) 

with X defined above. 

Finally, we calculate diagrams such as the one shown in Figure Th ese are 
vertex corrections to the "ordinary" second-order diagrams. Figure J. 7 shows 
the first part of Figure J. 6, up to the intermediate state C. The plan is to find an 
effective Hamiltonian for this half of the diagram (the actual vertex correction), 
and then use it to calculate the full diagram. Calling the new intermediate NN 
state _B, the first of the diagrams pictured in Figure J. 7 is 



B,Pb =■= 



1 



fioji + Ed, - Eb - P^/2md - E^+ie 



B,Pb\^ (ai ■ A{x,)) ^T^^e'-^--- \d,,P,,-i,). (J.77) 




Figure 3.7: Possible time-orderings for the first half of the meson-exchange 
diagram pictured in Figure p.6| (a). All are equal. 
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Next, we make the assumption of equation (J. 9) to simplify the denominator, 
and then remove the complete set of B states to get 



^f? = ^^^^{C,Pc\{~T2+n^+T2-n+)x (J.78) 

(fiq giq-{x2-Si) . 



TTeh-'f 



3 r2 



■ q) — — -2-7 i^i ■ ^(^1) \di,Pi,-fi ). 

m4 -I- fin^ V / 



(2^)3 ' ' ^' ml + h^q^ 
Performing the g-integration yields 

1 p f 2 ^ ^ / p-m„\x2~xi\/h 

T- = ^(fiPd[nx.4<^2-v( )x (J.79) 



where the gradient is with respect to a;2 — afi (which will become f). The term 
where the pion is destroyed, rather than created, at the photon vertex, can be 
calculated similarly, and turns out to be identical to equation ( J.79| ). Adding 
these two terms to the two terms where nucleons 1 and 2 are exchanged (see 
Figure J. 7 for the Feynman diagrams) gives 

,^i,tot en/2 ~i f/^ ^r^^\f^ e^m.|S2-Si|A 

= ^i^^Pc I [n X r4 I (a, . A(.0) [a. ■ V 

[<^2-A{x2)) U-V |^^_^^| )| M.,P.,7. >■ (J-80) 

Therefore, we can write the effective Hamiltonian which includes all terms where 
one of the photons interacts with a pion as 

H'" . _^p,xf.u|(,-,^^(ft))(<?.^V^5-^ 

(.^..A(f.))(.,.V^^-^)|. (J.81) 

This formula applies to either incoming or outgoing photons, since no assump- 
tions about the photons have been made. Since this Hamiltonian changes 
the deuteron isospin, we must use second-order perturbation theory with an- 
other isospin-changing Hamiltonian. The only operator available for this is t = 
(gi,-gn)/ 2. Th e largest contribution comes from i?'"* = - / A(-^\i)d^^ 
(equation G.16 ). Therefore, the vertex corrections (vc) to the scattering ampli- 
tude are 

W^' = (J.82) 



C,Pc 



huj, + Ed, - Ec~ P^/2md + ie 
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C,Pc 



-hujf + Ed^ - Ec- Pc/'^md + i 



■{C,Pc,l^,lf\H^'''\d^,P,,J,) 



fioji + Ed, - Ec- Pll^rtid + ie 
c,Pc 



( CFcl y J(^)(0 • \d^A,l^). (J.84) 



C,Pc 

1 



The calculations are very similar to the ones in Appendix G, so we give just 
one intermediate step in the calculation o f TJ^ '^^ . After inserting the expres- 
sions for the Hamiltonians (from equations J. 81 and G.16 ) and performing some 
simplifications, we get 



LcScJcMc W LML'M' ^ ^ 

+ i-lf^^n] D'i;?_,^ (0, -e, -<f>W2niL' + 1) ^^)x 

j rdrr'dr'j,,_,{^)u,{r')G,^{r,r';Eo)M^)j^^^ u,{r) x 

( /'11|| [YL'-i®t]L' WLcScJc )( I'llMfliidi ■ hM^2-f) + 

(-1)^(<?2 •r)}yLM laiM, ). (J.86) 

We have used the fact that 

( r = I [n X f2]. I T = 1 ) = -2i. (J.87) 



Equation (J.21) can be used to help the evaluate the second matrix element 
above. The final result is 

l+Jc l+L" L+L' oo oc 

^r = -EE E E E EE E E ^J-^^) 

1=0,2 I' =0,2 j=\i-j^ \ Jc = \l-L"\ L" = \L-L' \ L=0L' = 1 j=- 1,0, 1 J = - 1,0,1 



m„m. 



n(i') m a ^V- ^ ,o.~ , , ,^ / (i^' + l)(2L^' + l)(2Jc + l) 
^M'-A/O,-6',-0)(eAj~i(2J + 1)(2L + l)y — 
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1 



L' 



Jc 



1 1 






L 




,7 i 









o) 



Jc 



J 1 



I L L" \ f 1 L" J 

-j j J \ -j -i J \ -Mi -i i Mi 

Jc Jc 

II 1 }{ l'n\\ [Yl'-1 <8) t]L' WJcOJc )( Jell >L" ||/ ) X 

1 L" J 



DC /■•OO 



r r' dr' j,,_,{'^)u,{r')GUr,r'; Eo)jr,{^)^- w,(r). 



/o ^0 2 
The other correction terms are: 



2 dr r 



l + Jr 1+L" L+L' oc oc 

EEE E EEEEE -(J-^^) 

;=0,2i'=0,2 j=|i_Jc| Jc = |l-i"| L" = |I--i:-'| L=0L'=1 i=-l,0,l j=-l,0,l 



72ttc- J'-li ' XiuJi 



rripml 



l^L+L' ^_^-^L'+i+j-L"-Mf+Jc-Mi-\i 



[Mp + (-1)V] 



Yl,M,{kf){il^U{2J+lW{L + 1){2L' + 1){2L" + l)(2Jc + 1) 

1 L' L" \ / 1 L' L" 

J \-j Mf - Mi - Xi-i Mi + Xi + i + j-Mf ' ^ 

1 1 1 \ / 1 L" J 
j i Mf-Mi-Xi-i-j Xi + Mi-Mf 

Jc L 1 \ f 1 J Jc 

-M, - A, Xi Mi J y -Mf Mf-Mi- Xi Mi + Xi 

1 r 1 

Jc Jc }{ JcOJcW [Yl-1 ® t]L win ){ I'W Yl" \\Jc ) X 

1 L" J 



OO pOO 



^0 



dr r' dr' j,_,(^^^)u,(r)G,,(r,r';£;o)j.K^)^^^^«.'(r-'), 



1 + ,/c 1+L" L+L' oc OO 

^r = -EE E E E EE E E p-^^) 

;=0,2 ;'=0,2 ,7=|i_,7f,| Jc = |l-L"| L" = |L-L'| L=0L' = 1 i=-l,0,l j = -l,0,l 



rripml 



IJ.p + (-l)^'Mn 



I? 



M'lA,(O,-^?,-0)(eAj-.(2J+ 1)(2L + 1) 



(L' + l)(2i" + l)(2Jc + l) 



47r 



1 L L" 



1 \ / 1 1 1 



1 L" J 



j j J \ i+ j —j —Mf i Mf — i 





L 




)(i 





o") 


Jc 






Mf- 





X 
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1 /' 1 

Jc Jc }{ JcOJcW [Yl'-i O t]L' win )( ril Yl" II Jc ) X 

1 L" J ^ 

^ f°° ujtr' LU-r d e"'""'"/^ 



l+Jc 1+L" L+L' oo oo 

A^r = -EE E E E EE E E x(J-9i) 

i=0,2i'=0,2 J=|i_Jc| Jc = |l-i"| i"=|I--i:-'| i=0 L'=li=-l,0,l j=-l,0,l 

^z:'M'(^/)(eA,)-.(2i+ + l)(2i' + l)(2i" + l)(2Jc + 1) x 

1 L' L"\fl L' L" 

J\-j Mf - Mi- Xi-i Mi + Xi+i + j - Mf 

1 1 1 \ / 1 i" J 
j i + j Mf - Mi- Xi-i- j X, + M,-Mf 

1 L Jc \f Jc J 1 

-Mf Xi Mf-Xi )\\- Mf Mf-Mi- Xi Mi 

Jc Jc 

1 r 1 W Z'llll [F^.i ® t]L \\JcOJc ){ JcW Yl" \\l ) X 



X 



X 



X 



1 L" J 



rdrr'dr'u_,{^)u,{r')G,Jry;E',)M^)-^^ 



Appendix K 

Useful Formulas 



Tensor product definition: 



(K.l) 

Wigner-Eckart theorem: 

(K.2) 

Formulas for uncoupling tensor products: 

(Jij2j|| [Tk^®n,]k \\3i32j')= (K.3) 

V(2i + l)(2i' + l)(2fc + l)<| ii' i'' / ^(ii||TfeJ|ji')(i2||Tfe, 

fci ^2 A; 



(jII [^fe^^TfeJ, ||/)= (K.4) 

(_l)feWV2fcTlX:{ y ^/ J j(il|r,J|j)(j||T,j|/). 
Useful reduced matrix elements: 



(niniioH-i)-V '""'"",r"""" (o o)- (-) 

( ||i:.SV/;)= (K.6) 

l{2Lf + l)(2Li + l)(2Ji + 1)(2J/ + 1)(2L + 1){2L' + 1)(25/ + 1) 
V 47r 
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L Q L 



Lf L Li 





L, I ^ 

/3(2L/ + l)(2ii + l)(2Ji + 1)(2J/ + l)(2i + l)(2i' + 1) 



(K.7) 



(-1)^ 



27r 



^/ "S*/ ^/ 

Si Ji 

L' 1 i 



o' J 



{LfSfJfW [YL>®t]^ \\LM 



(K.8) 



/3(2i/ + l)(2^i + l)(2'/i + 1)(2J/ + l)(2i + 1)(2L' + 1) 



1677 



f Sj Jj 

[ L' 1 L 

{LfSfJfW ||ii5iJi) = 



L/ L' Li 




(K.9) 



Li{Li + l){2Li + 1)2(2% + l)(2Ji + 1)(2J/ + l)(2i + 1)(2L' + 1) 



47r 



i' Li 




L' 1 L 
Lj Lf Li 



Addition of two spherical harmonics: 



Y Y-, / = 



E(-i) 



JM 



M./ (2j + l)(2/ + l)(2J+l) /j / J 
47r V 



(K.IO) 

m m' -M '^-^^ 



Spherical Bessel recursion relations: 
2[+l 

{2l + l)j:{x) = 
The gradient formula: 

V$(r)riM(r) = 



+ ji+i{x), 
lj,-i{x) - {l + l)j,+i{x). 



(K.ll) 
(K.12) 



V X 



f{r)TLLM{f] 



2L + 1 \dr 

_d _ L 
rfr r 



+ ) Hr)TLL-lM{f), 



m^^^TLL+iMif) +(K.14) 



L + 1 



dr 



fir) 



L + 1 
2L+1 



Tl 



L-IM 



APPENDIX K. USEFUL FORMULAS 



158 



V X 

V X 

V- 
V • 



fir)TLL+lMif) 

f{r)TLL-lM{f) 
f{r)fLLM{f) 

f{r)fLL+iM{r) 

f{r)fLL_lMif) 



d L + 2 
dr r 

d L 



fir) 



L 



2L + 1 



TLLMif)lK.15) 



L+1 



= 0, 



(K.17) 



L + 1 



2L + 1 \dr 



d ^ L + 2 



L 



2L+1 \dr 
Vector spherical harmonic definition: 



d L-1 
-r + 



fir)YLMir), (K.18) 
fir)YLMif)- (K.19) 



f^.Mif)^ E (-1)— V27Ti(_^^ _i) 

!^=-l,0,l ^ ^ 



f-„YL^M+i,{f)- 

(K.20) 



